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INTRODUCTION 


The successful parametric modeling of the aerodynamics for an airplane operating 
at high angles of attack or sideslip is performed in two phases. First, the aerody- 
namic model structure must be determined and second the associated aerodynamic param- 
eters (stability and control derivatives) must be estimated for that model. Though 
the aerodynamic model structure is known to be linear at low angles of attack, the 
appearance of nonlinearities at higher angles of attack has been a prominent feature 
in several recent reports of flight test results (refs. 1, 2, 3). Since a large 
number of possible nonlinear terms could contribute to the aerodynamic function, some 
method must be developed that examines only the influential terms while ignoring 
those that are superfluous. One possibility is to look at all combinations of linear 
and nonlinear terms. However, as pointed out in reference 4, the number of models to 
be considered grows too fast with the number of possible terms for such a technique 
to be practical. The use of stepwise regression was suggested in reference 5. 

Stepwise regression examines each term as to its usefulness in improving the model 
(by reducing residual variance). Candidate model terms are added one at a time 
and/or deleted one or more at a time until no more candidate terms can pass a given 
test of statistical significance. This provides a least squares equation error 
estimate of the model structure and associated parameters at each step of the model 
building* 

The stepwise regression technique has been used (refs. 1, 2, 3) for analyzing 
flight data from high angle-of-attack and large amplitude maneuvers. The purpose of 
this paper is to document two versions of a stepwise regression computer program 
which were developed for the determination of airplane aerodynamic model structure 
and to provide two examples of their use. It is assumed that the reader is familiar 
with the airplane equations of motion. One should read references 1 and 2 for appli- 
cations of the technique to actual flight data. 

The two computer programs that are the subject of this report, STEP and STEPSPL, 
are written in FORTRAN IV (ANSI 1966) compatible with a CDC FTN4 compiler. Both 
programs are adaptations of a standard forward stepwise regression algorithm (ref. 6). 
The purpose of the adaptation is to facilitate the selection of an adequate mathe- 
matical model of the aerodynamic force and moment coefficients of an airplane from 
flight test data. The major difference between STEP and STEPSPL is in the basis for 
the model (found in SUBROUTINE DATASET). The basis for models in STEP is the 
standard polynomial Taylor’s series expansion of the aerodynamic function about some 
steady-state trim condition (see refs. 1 and 3). Program STEPSPL utilizes a set of 
spline basis functions (refs. 3 and 7). 

The paper is organized as follows. After this introduction is a section 
describing the approach and rationale of the program. The main program and sub- 
routines are each described as to their respective purposes and dimensioning informa- 
tion. Next, a section addresses the interpretation of output based on two examples. 
There are seven appendices. Appendix A is the listing for STEP. Appendix B is a 
listing of the NAMELIST/INPUT/f or the first example. Appendix C consists of the 
output for the first example (which demonstrates STEP). Appendix D is a sample job 
control deck for running the example in a batch mode at the Langley Research Center 
computer center. Appendix E is the STEPSPL listing. Appendix F contains the output 
of example 2, (which demonstrates STEPSPL). Finally, appendix G contains a sampling 



of options for the spline model basis used in conjunction with STEPSPL. The 
interested reader can start by running the given test case and then modifying the 
program to fit his specific use. 


STEPWISE REGRESSION 

This section describes the basic principles and features of the stepwise regres- 
sion which is used to determine aerodynamic model structure from flight data. It is 
assumed that the general structural form of the mathematical model for the aero- 
dynamic force and moment coefficients can be written as 


y ( t) = 0 Q + 9 1 x 1 (t) + 0 2 x 2 (t) + ... 9 n x n (t) 


(1) 


where 

y ( t ) aerodynamic force or moment coefficient ( c x' C Y' C Z' C m' C £' C n^ at time t 

x • ( t ) airplane state plus control variables ( a,q, B,p, r , 6 e , S a , <5 r ) and their com- 

binations at time t ( j = 1 , 2, . . . , n) 

0^ airplane stability and control coefficients ( j = 1 , 2, ..., n) 

0 Q constant reflecting and initial steady-state condition. 

The forward stepwise regression described in this paper begins with the assump- 
tion that there are no variables in the postulated regression equation other than the 
bias term 0 Q . An effort is then made to find an optimal subset of variables by 
inserting independent variables into the model one at a time. The first independent 
variable selected for entry into the equation is the one that has the largest cor- 
relation with the dependent variable y. Suppose that this variable is x 1 . This is 
also the variable that produces the largest value of the F-statistic for testing the 
significance of regression. The variable is then entered if the partial F-statistic 
of its associated parameter, 0^ exceeds a preselected critical F-value. 


F P = 


1 

2 A 


> F 


crit 


a 2 A 

where 0 is the estimated parameter associated with x 1 and s (0.) is the 
1 1 l 

variance estimate of 0^ 

The second variable chosen for entry is the one that now has the largest 
correlation with y after adjusting for the effect on y of x 1 . These correla- 
tions are referred to as partial correlations. In general, at each step, the in- 
dependent variable having the highest partial correlation with y is added to the 
model if the partial F-statistic of its associated parameter exceeds the preselected 
F - t . At each step of the procedure, all variables entered into the model pre- 
viously are reassessed by examining their corresponding partial F-statistics . A 
variable added at an earlier step may be redundant because the relationship between 
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it and the remaining variables now in the equation has reduced its value of F^ to 
less than F crit _. happens, the insignificant variable is deleted from the 

regression model. The procedure terminates when all significant terms have been 
included in the model. 


Five Associated Information Criteria 

At each stage of the stepwise regression, as a new variable enters the model, 
five useful quantities are calculated. All these quantities should be examined for 
the final model selection. First, the user can consider the total F-value for a 
given model of Q variables calculated as the ratio of the mean square due to the 
regression to the mean square of the residual. This ratio is given as 


l (Jin - y ) 2 

p = i=±_ _ n_-_J2 

N A Q - 1 

l (yd) - y(i)) 2 

i=1 


where 


- i N 

y = n E y(i) 

i=1 


This number usxially increases to some maximum value as new variables enter the 
regression, but then decreases slightly as the new terms are less effective in reduc- 
ing the residuals. Heuristically , the maximum F-value represents a model which best 
fits the data with a minimum number of parameters. Second, the squared multiple 
correlation coefficient R 2 is calculated. This number, expressed as a percentage, 
is a measure of the usefulness of the terms, other than 0 q/ in the model. The value 
of R 2 would be 100 percent for a model that perfectly fit the data. Third, at each 
stage, the partial F- values F D for each parameter are printed. The user should 
look for consistency in the value of F . For example, if one value of F p is only 
slightly greater than F rr * t and all other values of .F are much greater, the user 
may not want to include the variable with the small value of F^ in the model. The 
fourth aid in model selection is the estimated normalized autocorrelation function 
for the residuals. The estimate of the autocorrelation function at lag h is 
given by 


1 N ~ h 

W(h ) = — ) v(.i) v(i -f h) (h 0, 1, M) 

N - h . L 

1=1 


where h is the lag number and M is the maximum lag number, which is usually 
10 percent of N. For data sampled each At second, the time separation associated 
with lag h is h*At. The normalized autocorrelation function is calculated as 
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W(h)/W(0). This function should approach that for white noise with a value of 1 at 
zero lag and values of 0 at lags of 1 to M. In applications, when the value of F^ 
for a parameter makes the utility of an independent variable questionable, the con- 
tribution of that variable to the actual model structure can be assessed by observing 
the effect of the variable on the autocorrelation function of residuals. The fifth 
number that is useful is the standard error in the residuals, a, which is printed at 
each stage of the regression. 

One learns from experience that not all of the five criteria listed above are 
’’optimally" satisfied for any single model. However, the stepwise regression and its 
associated information criteria do significantly reduce the number of possible models 
from which the user must choose. Moreover, as the model structure is determined, so 
are the parameter estimates. Finally, ambiguity in the model selection can also be 
resolved by requiring that the estimated parameters make sense physically and that 
the selected model have good prediction capability. 


Selection of Candidate Model Variables 

The selection of a set of candidate model variables from which the stepwise 
regression can build a model should rely on the user's a priori knowledge of the 
physical system that is to be modeled. For the airplane, assumptions as to the most 
influential variables and symmetry considerations have led to the following logic for 
selection of candidate model variables for a spline analysis of the longitudinal 
maneuver. The range of the independent variable which is most important in the 
determination of the dependent variable is partitioned into several subsets, each 
having support on less of the range than the previous subset. For example, the force 
coefficient is mainly dependent on a. Hence if a = {z|a < z < b} , then the a 

range, [a,b] , is divided according to the spline basis functions as follows: 


(a - a. ) 

i 


m 

+ 


(a 


a. ) 
i 


m 


(a > a. ) 

i 


(a < a. ) 

i 


The values of are called knots. An example of the function is given in 

figure 1. The four knots in this figure are at a = 2°, 4°, 6°, and 8°. Hence, 

(a - a 1 )° = 1 for a > a 1 = 2°, and (a - a 1 ) ° = 0 
j g* i * + 

for a < . Similarly, 

(a - a 2^ ~ 1 for a * 4° , and (a - “ 0 ^ or a < 4° , 


and so forth, for the rest of the functions. If the order of the function, 

denoted by the superscript m is other than zero, say 2, then 


( a - a 1 ) 


2 

+ 


(a - a 1 ) 


2 

+ 


for a > a-j 


and 
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Though it appears lengthy and awkward, the following formulation of the FORTRAN 
code allows for simple deletion, addition, and/or change in candidate model vari- 
ables . 

DO 910 1=1 , NPTS 
X(1 , I ) =ALPH( I ) 

X(2,I)=C/(2*VEL(I) )*Q(I) 

X ( 3 , I ) =DELE ( I ) 

DO 911 111=4,39 
911 X(III,I)=0. 

IF (ALPH(l) .GE.XKN0T( 1 ) ) X( 4 , I ) =ALPH( I ) -XKN0T( 1 ) 

IF (ALPH(I) .GE.XKN0T( 1 ) ) X ( 5 , I ) =X ( 2 , I ) 

IF (ALPH(I) .GE.XKN0T(2) ) X( 6 , I ) =ALPH( I ) -XKN0T( 2 ) 

IF ( ALPH ( I ) . GE . XKNOT ( 2 ) ) X( 7 , I ) =X( 2 , I ) 

IF ( ALPH ( I ) . GE . XKNOT ( 3 ) ) X( 8 , I ) =ALPH( I ) -XKNOT ( 3 ) 

IF ( ALPH ( I ) . GE . XKNOT ( 3 ) ) X ( 9 , I ) =X ( 2 , I ) 

IF (ALPH(I) .GE.XKN0T(4) ) X( 1 0 , I ) =ALPH( I ) -XKNOT ( 4 ) 

IF (ALPH(I) .GE.XKNOT(4) ) X ( 1 1 , I ) =X (2,1) 

IF ( ALPH( I ) .GE. XKNOT ( 5 ) ) X( 1 2 , 1 )=ALPH( I ) -XKNOT ( 5 ) 

IF ( ALPH ( I ) . GE . XKNOT ( 5 ) ) X ( 1 3 , I ) =X ( 2 , I ) 

IF ( ALPH ( I ) . GE . XKNOT ( 6 ) ) X( 1 4 , I ) =ALPH( I ) -XKN0T( 6 ) 

IF ( ALPH ( I ) .GE . XKNOT ( 6 ) ) X ( 1 5 , I )=X( 2 , I) 

IF ( ALPH ( I ) . GE . XKNOT ( 7 ) ) X( 1 6 , I ) =ALPH( I ) -XKN0T( 7 ) 

IF ( ALPH ( I ) . GE . XKNOT ( 7 ) ) X ( 1 7 , I ) =X (2,1) 

IF (ALPH(I) „GE.XKN0T(8) ) X( 1 8 , I ) =ALPH ( I ) -XKNOT( 8 ) 

IF ( ALPH ( I ) . GE . XKNOT ( 8 ) ) X( 1 9, I )=X( 2,1) 

IF ( ALPH( I ) .GE. XKNOT ( 9 ) ) X( 20 , 1 ) =ALPH( I ) -XKNOT ( 9 ) 

IF ( ALPH ( I ) . GE . XKNOT ( 9 ) ) X ( 21 , I ) =X( 2 , I ) 

IF (ALPH(I) .GE.XKNOT(IO) ) X( 22, I ) =ALPH( I ) -XKNOT ( 1 0 ) 

IF ( ALPH ( I ) . GE . XKNOT (10)) X( 23, I )=X( 2 , I ) 

IF (AIiPH(l) .GE. XKNOT ( 11)) X( 24, 1 ) =ALPH( I ) -XKN0T( 1 1 ) 

IF (ALPH(I) ,GE.XKNOT( 1 1 ) ) X ( 25 , I ) =X ( 2 , I ) 

IF (ALPH(I) .GE.XKNOT(1 2) ) X( 26 , I ) =ALPH( I ) -XKNOT( 1 2 ) 

IF ( ALPH ( I ) . GE . XKNOT (12)) X( 27 , I ) =X( 2 , I ) 

IF (ALPH(I) .GE.XKNOT(13) ) X( 28 , I ) =ALPH( I ) -XKNOT ( 1 3 ) 

IF ( ALPH ( I ) . GE . XKNOT (13)) X( 29, I )=X( 2 , I) 

IF (ALPH(I) ,GE.XKNOT(14) ) X( 30 , I ) =ALPH( I ) -XKNOT ( 14 ) 

IF ( ALPH( I ) .GE. XKNOT ( 14)) X( 31 , I ) =X ( 2 , I ) 

IF ( ALPH ( I ) . GE . XKNOT (15)) X( 32 , I ) =ALPH( I ) -XKNOT ( 1 5 ) 
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910 CONTINUE 


IP ( ALPH ( I ) . GE • XKNOT (15)) X( 33 , I ) =X( 2 , I ) 

IF ( ALPH ( I ) * GE • XKNOT (16)) X ( 34 f I ) =ALPH( I ) -XKNOT( 1 6 ) 
IP ( ALPH ( I ) • GE • XKNOT (16)) X( 35 , 1 ) =X( 2,1) 

IF ( ALPH ( I ) • GE • XKNOT (17)) X ( 36 , I ) =ALPH ( I ) -XKNOT (17) 
IF ( ALPH ( I ) . GE • XKNOT (17)) X( 37 , I ) =X ( 2 , I ) 

IF ( ALPH ( I ) . GE • XKNOT ( 7 ) ) X ( 38 , I ) =X ( 3 , I ) 

IF ( ALPH ( I ) • GE • XKNOT (13)) X( 39 , I ) =X ( 3 , I ) 


In the preceding printout, VEL(l) = airspeed V at t^ , Q = pitch rate q, 
NPTS = number of data points N, C = wing mean aerodynamic chord c , and 
X(J,I) = value of jth model variable at t^. The symbols XKNOT( ) indicate knots 
Cor specific values of a. The code above actually gives the logic for creating the 
(39 x n) matrix containing the time histories of each of the 39 candidate independent 
variables. The 17 knots in angle of attack can be set at any value the user deems 
adequate for the data by setting XKNOT(I) in the program with I = 1,17. Changing 
the candidate model variables can easily be accomplished by substituting the new 
variable for any of the 39 candidate variables listed. The number of candidate vari- 
ables is limited only by the size of the computer memory. 


DESCRIPTION OF PROGRAM STEP 
Main Program 

The main program for STEP is dimensioned to accept measurements of flight data 
at 500 time points. For example, at a data rate of 20 points per second, this dimen- 
sioning allows 25.0 seconds of data to be used. The main program includes all logic 
for the actual regression procedure as well as most of the printing logic. The name- 
list INPUT is read in the main program. This namelist contains data for airplane 
mass, geometry and inertia characteristics, initial conditions, option switches and 
starting time. The namelist will be discussed in detail below. The main program 
also does all calculations involving the correlation matrix and analysis of variance. 
It provides for the printing of the partial F-values for variables in the regression, 
the estimates of their coefficients and standard errors of those coefficients. The 
total F- value for the model, the percent variation from the mean explained by the 

regression model, and the variance of the residual sequence are also printed from the 

main program. 

Dimensio ns - As mentioned, most data arrays are dimensioned 500 to accommodate a 
maximum of 25 seconds of data at a rate of 20 points per second. However, this and 
other dimensions may be adjusted by the user to conform to individual computer 
capability. This section is written to aid the user in such changes. Let MAXNPTS be 
the maximum number of data points to be analyzed and N-1 be the maximum number of 
independent variables to be considered for the regression. For example N-1 is 24 
in STEP since the lateral equations have 24 candidate model variables as seen from 

the X array in lines 79 through 102 of subroutine DATASET (appendix A). The dimen- 

sions for arrays in STEP are as follows: 

T, Y, YHAT , XNU , AX, AY, AZ , PDOT, QDOT, RDOT, VEL, P, Q, R, and 

QQ are each dimensioned MAXNPTS 

X is the two-dimensional data array and should be dimensioned N x MAXNPTS 

S and XXSUM are two-dimensional work arrays and are dimensioned N x N 
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ICNT, IORD, B, V, STDER, FPART, TN, PPLT, FPLT, XSUM, XBAR, SIGMA, 
PPRT are each dimensioned N 


W and XLAG are arrays for the autocorrelation function and its lag number* 
They should be dimensioned at least * MAXNPTS/1 0 . 

APR is used in the calculation of the PRESS (Prediction Error Sum of Squares) 
criterion* It is a two-dimensional array with dimensions N x ( MAXNPTS + N) . 

RR, A, AP are two-dimensional arrays containing variance and covariance 
information and are dimensioned (2N-1) x (2N-1). 

Namelist Input - A namelist called INPUT communicates airplane geometry, mass 
and inertia characteristics as well as initial conditions and logic switches for 
program options. The elements of INPUT and their definitions are listed alpha- 
betically as follows: 


ALPHT - angle of attack trim value (radians) 
BETT - angle of sideslip trim value (radians) 
BSP AN - wing span (meters) 

CBAR - wing mean aerodynamic chord (meters) 


DELAT - aileron displacement at trim (radians) 

DELET - elevator deflection at trim (radians) 

DELRT - rudder displacement at trim (radians) 

FCRT - critical F-value for entry or elimination of a term in the model 

A nominal value between 5 and 10 for FCRT has proven effective from experience 
with high angle-of-attack airplane data. 

G - acceleration due to gravity (m/sec 2 ) 


IACEL0P - option switch to read angular accelerations from data 

0 - read from data 

1 - calculate by cubic spline subroutine 


IEQN - 


indicates which equations are 
If LATOP = 0 then IEQN = 1 

2 

3 

If LATOP = 1 then IEQN = 1 

2 

3 


to be fit: 
for C 
for C z 
f °r C m 
for C 
for C p 
for r 


IFILOP - option switch to incorporate low pass filter on lateral acceleration 
measurements 

IFILOP = 1 for filter active 
IFILOP = 0 for filter inactive 
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I FLAG 


- option switch to have first LINMAX terms considered before any other 
terms. LINMAX is set in the main program. 

I FLAG = 1 activates the option 

I FLAG = 0 all terms are searched from the first pass on 

LINMAX is set to 3 in line 38 of Program STEP (appendix A) for the longitudinal 
option (LATOP = 0) and 5 in line 40 for the lateral option (LATOP = 1). These values 
allow for consideration of the first three candidate model variables a, q, and 6e 
for the longitudinal equations and the first five candidate model variables $, p, 
r, 6a, and 6r for the lateral equations. To consider any J terms first, in 
general, Linmax should be set to J and those J terms should be the first J 
terms entered into the two-dimensional X array in SUBROUTINE DATASET. 

I PLOT - option to activate plotting 

IPLOT = 1 activates plotting 
IPLOT = 0 for no plotting 

IPRESOP - option to invoke PRESS calculation 
IPRESOP = 1 activates option 
IPRESOP = 0 for no PRESS calculation 

IPSKP - For IPRESOP = 1, IPSKP selects every (IPSKP)-th point for calculation 
of PRESS 

ITRIMOP - option switch to read trim values from first data point to be 

analyzed. If ITRIMOP = 1, the namelist supplied values (or default 
values of 0) for ALPHT, BETT, DELAT, DELET, DELRT are used. 

IF ITRIMOP = 0, the values of angle of attack, angle of sideslip, 
aileron deflection, elevator deflection and rudder deflection at time 
TS are used for ALPHT, BETT, DELAT, DELET, DELRT, respectively. 

2 

IX - moment of inertia about longitudinal body axis (kg m ) 

IXZ - product of inertias (kg r) 

2 

IY - moment of inertia about lateral body axis (kg m ) 

2 

IZ - moment of inertia about vertical body axis (kg m ) 

LATOP - option switch for lateral equations 

LATOP = 1 for fitting lateral equations 
LATOP = 0 for fitting longitudinal equations 

M - airplane mass (kg) 

NEQ - the number of equations to be fit (as opposed to IEQN which indicates 
which NEQ equations are to be fit). NEQ can be 1, 2, or 3. 

NPTS - number of data points to be fit 

PT - roll rate trim value (rad/sec) 

QT - pitch rate trim value (rad/sec) 

o 

RHO - atmospheric density (kg/m ) 
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RT - yaw rate trim value (rad/sec) 
SAREA - wing area (m 2 ) 

TS - starting time for data to be fit 


Subroutines 

SUBROUTINE DATASET - The purpose of SUBROUTINE DATASET is to read flight data 
from a file (TAPE 1 ) into the program and to set up the data array of time histories 
of the candidate independent model variables . The candidate variables for both the 
longitudinal and lateral options are given in Table 1 . The user can easily change 
any of the candidate model variables in the table to meet his own needs . The 
candidate model variables presented here have proven to be useful in the work 
reported in reference 1 . 

The number of candidate model variables is limited only by the size of computer 
memory. Any changes in the number of candidate model variables in the X array of 
SUBROUTINE DATASET should be reflected in the value of N in the main program and 
possibly the dimensions which depend on N throughout the program. 

SUBROUTINE DATASET also calls a cubic spline differentiation subroutine 
(SUBROUTINE SECDER AND FUNCTION DERSP ) for the calculation of angular accelerations 
from the measured angular rates. These calls can be eliminated as can the SUBROUTINE 
SECDER and FUNCTION DERSP if the user has measured angular accelerations available. 

SUBROUTINE AUTO - The normalized autocorrelation function for the residual 
sequence is calculated. XLAG and W must be dimensioned at least MAXNPTS/10. 

SUBROUTINE FIL - This subroutine is a low pass filter for the smoothing of data 
in the time domain. The algorithm is taken from reference 8. When this filter 
routine is active (IFILOP = 1) the user must choose FC and FT (which define the 
frequency range in Hz for band pass roll off) in this subroutine so that H(I), 

1=1, NPTS/2 is always defined. 

Subroutines for PRESS calculation - STEP and STEPSPL programs use six sub- 
routines for PRESS calculations . The main program calls upon three primary routines: 
PRESS, UPDATE, and PSET. These in turn call on three secondary routines: REDEF , 

INTRCHG, and RANDOM. Subroutines PRESS and UPDATE are called during each pass as 
model variables are added or deleted. The PRESS routine simply computes the value of 
PRESS associated with each candidate variable. This is done without any effect on 
the regular stepwise regression calculations. Subroutine UPDATE is used to modify 
the normal equations to reflect the change in model variables during each pass. A 
separate set of normal equations is used by PRESS so that the stepwise regression and 
PRESS calculations can proceed independently. 

Subroutine PSET is called once at the start of a run to establish the dataset to 

be used in the PRESS computations . As described in reference 1 , for a large number 

of data points PRESS approaches the residual sum of squares (RSS). Therefore it may 

be necessary when handling large datasets (number of points greater than 100) to use 
a reduced number of data points. The IPSKP variable controls the number of data 
points to be used by PRESS. The selection of 30 to 40 points has proven to be best. 
IRAN =1 is the flag which indicates the reduced dataset is to be randomly selected. 
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The secondary routines RANDOM, REDEF , and INTRCHG provide a few simple opera- 
tions. RANDOM returns a uniformly distributed random variable which is used to 
randomly select data when required by PSET. Subroutine INTRCHG is used by UPDATE to 
interchange rows and columns in the normal equation matrices* Subroutine REDEF is an 
initializing routine used to prepare the appropriate matrices for PRESS computations. 


DESCRIPTION OF PROGRAM STEPSPL 

STEPSPL differs from STEP mainly in the dimensions of arrays and in the 
SUBROUTINE DATASET. Program STEPSPL is dimensioned to accept data lengths of 900 
points (45 seconds of data at 20 points per second). There is provision for 39 
independent variables which will be spline "+" functions and for 23 spline knots. 

The " + " function is defined as (Aa)™ = (a - = (a - a ^) + if a * (Aa) + “ 

if a < a k , where a k is the value of the kth knot in angle of attack in radians. 

Example 2 demonstrates the use of STEPSPL. Appendix E contains the STEPSPL 
listing. 

Dimensions - Let MAXNPTS be the maximum number of points to be analyzed and 
let n — 1 be the maximum number of independent variables to be considered. The 
dimensions of arrays are as follows: T, Y, YHAT , XNU, AX, AY, AZ, PDOT, 

QDOT, RDOT , VEL, P, Q, R are each dimensioned MAXNPTS. 

X is the two-dimensional data array and is dimensioned N x MAXNPTS. 

S and XXSUM are two-dimensional work arrays and are dimensioned N x N. 

ICNT, B, V, STDER, FP ART , XSUM, XBAR, and SIGMA are each dimensioned N. 

W and XLAG are arrays for the autocorrelation function and its lag number. 

They should each be dimensioned at least MAXNPTS/10. 

RR, A, and AP are two-dimensional arrays containing variance and covariance 

information and are dimensioned (2N - 1) x ( 2N - 1 ) . 

Namelist Input - The namelist for STEPSPL is the same as that for STEP with the 
exception of the variables IPRESOP, and IPSKP which apply to PRESS subroutines not 
found in STEPSPL. 


USING STEP AND STEPSPL 
Aids in the Selection of an Adequate Model 

Since there is no cost function which ensures that the best model has been 
found, STEP and STEPSPL provide a subset of all possible models. From this subset, 
one must make the selection of an adequate model for the data at hand. To assist in 
the selection of an adequate model, the programs provide several statistical and 
informational parameters at each step of the fitting process. These parameters are 
as follows: 

1 . The partial F-values for the coefficients of all variables that are 
currently in the model. 
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2. The total F- value associated with the current model. 

3. The square of the correlation coefficient in percent corresponding to the 
percent variation from the mean of the data that is explained by the current 
model. 

4. The Prediction Sum of Squares (PRESS) criterion - The scalar PRESS , 
corresponding to the £th subset of model variables, is defined as 


A 2 
PRESS = l (y(i) - y[i/x(1), ..., x(i-1), x(i+1), x(N)] } 


where y(i) is the ith response of the system and y(i/...)^ is the least 
squares estimate of E{y(i)} for the £th subset. A Note the ith observa- 
tion, x(i), is not used in forming the estimate y(i/...)£. The model 
corresponding to the smallest value of PRESS is the best predictor model. 

It is also a parsimonious model since PRESS reflects the added cost of 
redundant model variables. 

5. The standard deviation of the residual. This should approach that 

calibrated for the instrument used to measure the dependent variable. 

Also at each point in the selection process the user is provided with a synopsis 
of variables currently in the model as well as the estimates of the coefficients of 
those variables and the standard error of those estimates. 


Example 1 

In this example STEP is run on a simulated data set. The program listing is 
found in appendix A. The simulated data set to which the program STEP is applied is 
a subset of the time history for the mathematical model given in figure 2. The 
subset consists of the 43 points corresponding to an angle of attack in the range 
14° < a < 16°. The true values for the parameters in this range are 



0.700 



-3.00 



-1 .00 


C = 0.0 
X 

q 


C = -20.0 
Z 

q 


C = +15.0 
m 

q 


C = 0.05 
X 6e 



- 1.10 



-1 .00 
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The namelist/INPUT/f or this example is found in appendix B and the output listing is 
in appendix C. In examining the namelist/INPUT/it is seen that: 

1. IPRESOP ~ 0 - No PRESS calculation will be made; with PRESS inactive, 
processing time is cut by about a factor of 2* 

2* TS = 0.0 - The first data point to be considered for fitting is that 
corresponding to time 0.0 sec on the data tape. 

3. NEQ = 3 - All 3 equations corresponding to the choice of LATOP will be fit. 

4. IEQN = 1, 2, 3 - The NEQ equations to be fit are 1, 2 , and 3. 

5. NPTS = 43 - 43 datapoints after TS are to be fit. 

6. I PLOT = 1 — The program will plot the measured and computed time histories, 
residual sequence, and autocorrelation sequence at each step of the 
regression. If PRESS is active, IPLOT = 1 will also allow for the plotting 
of a synopsis of F-values and PRESS values after the last significant variable 

has been added to the model. 

7. IFLAG = 1 — Selects the option whereby the first LINMAX terms are considered 
before any others. These terms correspond to a linear model. 

2 

8. SAREA = 13.74 - The wing area Cor the airplane is 13.74 m . 

9. BSPAN = 9.98 - The airplane wingspan is 9.98 m. 

10. CBAR = 1.40 - The airplane wing mean aerodynamic chord is 1.40 m. 

11. M = 1055 - The airplane has a mass of 1055 kg. 

12. RHO = 1.0272 - The mean atmospheric density during the maneuver was 
1 .0272 kg/m 3 . 

13 . g = 9.81 m/sec 2 - is the gravitational acceleration constant. 

14. ix = 2357 - The moment of inertia about the longitudinal body axis is 
2357 kg m 2 . 

2 

15 . xy = 3051 - The moment of inertia about the lateral body axis is 3051 kg m . 

2 

16. IT. = 4833 - The moment of inertia about the vertical body axis is 4833 kg m . 

2 

17. IXZ = 177. - The product of inertia is 177 kg m . 

18. DELET = -0.08318 - The elevator displacement at trim initial conditions is 
-0.08318 rad. 

19. ALPHT = 0.2095 - The trim angle of attack is 0.2095 rad. 

20. BETT = 0 - The trim angle of sideslip is 0. rad. 

2 i DELAT = 0 — The aileron displacement to trim is 0. rad. 
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22. DELRT = 0 - The rudder displacement to trim is 0. rad. 

23. QT = 0 - Trim pitch rate is 0. rad/sec. 

24. PT = 0 — Trim roll rate is 0. rad/sec. 

25. RT = 0 - Trim yaw rate is 0. rad/sec. 

26. FCRT = 5 - The critical partial F- value for entry into the regression is 5. 

27. ITRIMOP = 1 - The trim values provided in this namelist will be used as 
opposed to the values of the associated variables at time TS. 

28. IPSKP » 10. - Indicates the the PRESS option, if activated by IPRESOP = 1, 
should select every 10th point for the evaluation of the PRESS. 

29. LATOP = 0 - Indicates that the longitudinal equations are to be considered 
for the fitting. 

30. IACELOP = 0 - Indicates that angular accelerations will be read directly from 
the data string. 

31. IFILOP = 0 - Indicates that the low pass filter will be inactive. 

After the namelist, the trim values for angle of attack, angle of sideslip, 
aileron, elevator, and rudder are printed. The output is continued in appendix C 
with a line of header information and a run identifier from the data tape, the value 
of IEQN and the number of points to be fit printed. It is seen that for RUN 1, equa- 
tion (1) (C x since the longitudinal option is active), is to be fit for 43 data 
points (NPTS =43). If an even number of points is specified for NPTS in the name- 
list, that number will be decreased by 1 by the program so that NPTS is always odd. 

Next is a listing of the relevant data for the points being fit. Here, for the 

longitudinal option, time, velocity, angle of attack, pitch rate, and elevator 
deflection are listed. If the lateral option had been chosen (LATOP =1), then time, 
velocity, angle of sideslip, roll rate, yaw rate, aileron deflection, and rudder 
deflection would have been printed. 

The next line indicates that the highest correlation between the measured depen- 
dent variable (a x here since LATOP = 0 and IEQN = 1) and an independent variable 

is for the first model variable, which is angle of attack. The partial F-value for 

this variable is 878, and its entry accounts for 95 percent of the variation. The 
standard deviation of the residual sequence (a - a ) is 

Measured x computed from model 
0.00139. The total F-value for this model is 857. 

The next line gives the least squares estimates for parameters currently in the 
model. The order of the estimates is the same as the entry of data into the X 
array in SUBROUTINE DATASET. Below the parameter estimates is found the estimated 

/v 

standard error for that estimate. Here, C v =0.651 with o„ =0.022 and the 

« C X n( 

model for C‘ x is at this stage: 
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C X " + C xJ a - a o> 


= 0.000229 + 0.651 (a - a Q ) 

An optional visual aid is provided in the form of a plot. An example of the plot and 
its interpretation is given for the STEPSPL program in example 2. 

The next variable chosen for the regression is variable 3, with a partial 

F- value of 5.59 x 1 0 8 . Since the entry of this term explains essentially 100 percent 
of the variation, it enhances the partial F-value of variable 1 also. The standard 
deviation of the residual sequence is now 0.38 x 10 8 and the F-value for the model 
is 5.98 x 10 9 . 

The new parameter estimates are — 0.700 and — 0.050. The respective 

A a fie 

standard errors are 0 =0.61 x 1 0 9 and O = 0.20 x 10 and the model is 

C V 


C v = C v + C Y (a - a ) + C (fie - 6e r o) 
x X o X a ° X fie 

= 0.000063 + 0.700 (a - a Q ) + 0.05 (fie - fie,o) 

This completes the fitting of the equation. 

Next, RUN 1, equation (2) is to be fit. Equation (2) (for LATOP =0) 
corresponds to the C z force coefficient. Again 43 points are to be fit. The most 
significant of the first three variables (since I FLAG = 1, the first LINMAX = 3 
are considered) is variable 2, q. With 58.76 percent of the variation explained, 
the model at this point is 


C Z “ C Z 0 + C Z q 2V 


Next variable 3, 6e, is added as being most highly correlated to the residual 
sequence of the previous model. With entry of the <$e term, 72.64 percent of the 
variation is explained and the model is now 


C Z - C Z Q + C Z q 2V T 


— + c„ (Se - Se,o) 


= -1.34 - 17.8 _ 0.705 (6e - 6e,o) 
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With the entry of the variable 1, a, 100 percent of the variation is explained and 
the final model is simply the complete linear model. 


C 7 = C„ + C 7 (a - a„) + C 7 + C*. (6e - 6e,o) 

o a 


T % 2V T ^ 


= -1.21 - 3.00 (a - a Q ) - 20.0 ~ - 1.10 (6e - Se,o) 


The third and final equation to be fit in example 1 is the pitching moment equation. 
The process goes as in the equations (1) and (2) with the first three terms incorpo- 
rated into the model. With the linear model completed, the model equation is 


= -0.731 - 1 .05a + 15.3 ^ - 0.996 ( Se - 6e,o) 


and explains 99.9 percent of the variation. However, the program adds, in the next 
step, variable 7. Note that the partial F-value of 29 for this variable is much less 
than and totally out of line with the first three terms (with partial F-values of 
2960., 5530., and 20,300.). Hence, the user might consider this last term to be 
superfluous and retain the linear model from the previous step. The unexplained # 
0.06 percent has been contributed by the spline differentiation of q to obtain q. 
When data for accelerations were read directly from the simulation program, 100 per- 
cent of the variation was accounted for. 


Example 2 

This example demonstrates the use of STEPSPL. STEPSPL is the basic STEP program 
with some dimension changes to allow for longer data lengths and the spline basis 
functions incorporated into SUBROUTINE DATASET. The simulated data for this example 
was generated by numerically integrating a model given by: 


C x = -0.180 + 0.700a + 0.050 ( 6e - 6e,o) 

C z =0.112 - 5.00a + 2.00 (a - 0.2269) + + 1.50 (a - 0.3142) + 

10.0 (a - 0.2269)^ „ 10. 0 (a - 0.3142)^ ~ 

- Se,o) 


- 10.0 


qc 

2V 


- 0.800 (5e 
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C m = 0.105 - 0.400a - 0.600 (a - 0.2269) + - 1.00 (a - 0.3142) + 
- 15.0 + io.o (a - 0.1745)° || + 10.0 (a - 0.2269)° || 


- 10.0 (a - 0.3142)° 


- 2.00 (fie - 5e,o) 


This longitudinal model is integrated using an elevator doublet input string* Random 
noise (with Gaussian distribution, zero mean, and a = 0.003) was added to pitch 
rate q. This a corresponds to ground calibration measurement error for the pitch 
rate gyro in previous flight testing. This a should yield a a of 0.08 for q 
and 0.023 for C m . 

Appendix E gives the results of applying STEPSPL to the noisy data sequence. 

The first information written is the namelist so that the user may quickly confirm 
that all elements on the list are correct. In this example, it is seen that only one 

equation (NEQ =1) is to be fit. That is the third equation ( IEQN = 3) which is the 

pitching moment equation since LATOP = 0. With IACELOP = 1, the spline subroutines 

are called to numerically differentiate the angular rate in order to derive angular 
accelerations. Since there is no PRESS option with STEPSPL, the PRESS associated 
options of example 1 are absent. Otherwise all INPUT options are the same as in 
example 1 . 

After the namelist, there is a listing of the angles of attack corresponding to 
cardinal knot positions. Following the knot values is a listing of trim conditions 
that the program will be using. Next is a line of header information giving a run 
identification number, the equation to be fit and the number of points to be fit. 
Here, it is seen that RUN 1, Equation 3 (C m ) is to be fit and 239 points will be 
used. Next the relevant data is printed. Since the longitudinal option has been 
selected, these data are time, airspeed, angle of attack, pitch rate and elevator 
deflection. 

Following the data listing is the actual fitting information. The overall list- 
ing of information is the same as described in example 1 for STEP. The major differ- 
ence is in the number of candidate model variables (which is now 39 plus a bias 
term). The maximum partial F-value (349.) is associated with the variable number 3, 
elevator deflection. The percent variation explained by the addition of this vari- 
able is 59.49 percent leaving a standard deviation of the residual sequence of 0.07. 
The parameter C m is estimated as -2.31 with a standard error of 0.12. The bias 
term, C , is estimated to be -0.000723. The next term selected by STEPSPL is vari- 

m Q 

able 1, angle of attack. The listing now shows that variables in positions 1 and 3 
are in the regression (by the H 1 1 s" in those positions). The percent variation 
explained by this model is 89.65 percent leaving a standard deviation of the residual 
sequence of 0.036. The total F-value for this model is 1022. The parameter 
estimates are: 
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0.246 



C = -1 .05 (±0.040) 
m a 

C m = -2.09 (±0.063) 
m Se 


The process continues until the entry of variable 13 corresponding to q 
(a - 0.1571)° the knot at a = 9° (9° = 0.1571 rad). The final estimated model is 
given in Table 2 and in figure 3, where it is compared with the true model used to 
generate the data. Though some of the values given in table 2 appear to be bad, it 
is seen from figure 3 that the overall model is very good. The program has approxi- 
mated the one knot in C at 10° by two knots: one is at 9° and the other is at 

m q 

11°. Thus, the table of values does not offer as good a feeling for the model as 
does the figure. 

In addition to the printout that has been discussed for this example and 
example 1, a plotting option is available. The subroutines used in the listings 
provided for STEP and STEPSPL are local to the LaRC computer center, but the user may 
combine whatever software is at his disposal to plot the same information. Figure 4 
contains the plot output for example 2. At each variable entry, three plots are 
generated. For example, figure 4(b) represents the entry of variable 3, <5e . The 

bottom plot in figure 4(b) displays the measured C m (t’s) and the C m computed by 
the model (solid line) at this stage of the regression. Above that, is a plot of the 
residual time history and the top plot is the autocorrelation function of the 
residual sequence. It is seen in figure 4(b), that the one variable model leaves 
quite a bit of structure in the residual sequence. In figure 4(c), the model, 
residuals, and autocorrelation sequence for the model containing Se and a are 
plotted. The autocorrelation sequence and the residual sequence are improved 
dramatically over figure 4(b). By figure 4(f), the visual aid displays a good fit, 
and good autocorrelation for the residual sequence. Figure 4(f) corresponds to the 
model containing 5 variables plus a bias term. The remaining parts of figure 4 all 
indicate a good fit and acceptable autocorrelation function. In general the plots 
have several applications to the curve fitting problem. One application is through 
structure that is left in the residual sequence. The user can look for new candidate 
model variables that might remove that structure. Secondly, if the structure is too 
fine for the user’s eye, the autocorrelation function may indicate that that struc- 
ture is present. Third, of course, is simply a picture of how well the computed 
curve fits the measured data. This also demonstrates the noisiness of the data. For 
example, a large variance in the residual sequence may indicate some filtering is 
required on the measured data. The indication for filtering is especially strong 
when the model fits the overall trends in the measured data but the residual variance 
is still large. 
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CONCLUDING REMARKS 


Two versions of a stepwise regression computer program which were developed for 
the determination of airplane model structure from flight data have been presented. 
The use of the program STEP with a Taylor's series expansion of the aerodynamic force 
and moment coefficients was demonstrated in example 1 . It is recommended that this 
program be used in regions where the variations in angles of attack and sideslip are 
not large but nonlinearity or aerodynamic coupling is suspected. Secondly, an 
example employing program STEPSPL was given. This program uses spline basis func- 
tions for the aerodynamic force and moment coefficients. It is recommended that 
STEPSPL be used when maneuvers having large variations in angle of attack and/or 
angle of sideslip need to be analyzed. The appendices contain the program listings 
and output for the two examples. 
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APPENDIX A 


This appendix contains a listing of PROGRAM STEP. 
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PROGRAM STEP 


74/74 OPT-1 


FTN 4.6*528 


82/10/28. 07.41.22 


. PROGRAM STEP (INPUT, OUTPUT, TAPE1,TAPE2, TAPES- INPUT, TAPE6-0UTRUT) 

REAL N,ir,IX,IZ,IXZ 

DIMENSION RR(49,4<?), A <40, 49), AP (49,49) 

OIMENSICN T(500 )> Y( 500 J, YHAT( 500), XNU(500) 

DIMENSION S(25, 251,3(25), V(25) 

OIMENSION ST0ER(25),FPART(25) 

DIMENSION W(50),XLAG(;0) 

DIMENSION I EON { 3 J ' 

DIMENSION TN(25J,PPLT(25)»FPLT(23J,TITLE(8) 

CCMMON/STAPT/ X ( 25 , 500 ) , XXSUH ( 25, 25 ) , XSUM ( 25 >, X8 AR { 23 ) , SIGMA <25 1 
COMMON /ACDATA/ S AR E A, R S P AN, C B AR, M, RHO, G, IX, I T, IZ, IXZ, DELET, ALPHT 
1»PETT*DELAT,DELRT,0T,PT,RT 
COMMON /AOP/ APR(25, 525), 00(500), PPRT ( 25 ) , PRSMIN 
COMMON /FLAGS/ IPSKP,NPTS,IDIM,JOIM,NMAX,IMIN,ICNT(25),IORD(25) 

* , IPPTS,LATOP,ITRIMOP,ICALL, iacelop,ifilop 
C fWTN/OPDER/ l£o,N 

COMMON/ ACC FL/ AX (500), A Y ( 500 > , A Z ( 500 > , P DOT ( 500 ), QOOT ( 500 ), ROOT ( 500 
ll»VEl(5CO),P(500)»0(5CO),R(500) 

EQUIVALENCE ( A » R R ) 

NAMELIST /INPUT/ IPRESOP,TS,NEO, IEQN, 

* NPTS, IPLOT, IFLAG, 

* SAPEA,BSPAN,C3AR, 

* M,RHO,G, 

* IX, IY, IZ, IXZ, 

* DELET, ALPHT, B ETT, DEL AT, DELRT, QT, PT, RT, 

* FCRT, ITRIMDP,IPSKP,LATO°,IAC£LOP,IFILOP . 

ALPHT»BETT»DELET«DELAT»0ELRT*OT*PT«RT*O. 

CALL PSEUDO 

l TOL •! . 0E-08 S I C AL L*0 ' 

READ(5, INPUT) •* - 

IF ( EOF (5 ) ) 2501,2503 ' - 

i VRITE(6, INPUT) ■ 

PRINT SUMMARY TITLES ON TAPE2 

WRITE (2,980) 

) FORMAT (1CX,*STEPVISE REGRESSION SUMMARY*) ' 

N-15 *LIN, MAX-3 

IF ( L AT OP. EC. 1 ) N"25 

IF(L ATOP. £0,1) L I N M A X » 5 


o ■« 
c > 
->■© 
r- m 


CO 1000 L* 1 , NEC 



PROGRAM STEP 


74/74 OPT-1 


FTN 4.84528 


82/10/28. 07.41.22 


IEO-IEON(L) JNGO-O 

ICALL-ICAU+l „ . .. .. ... 

CALL OATASETITS* T,Y,X) _ „ 

WRITE ( 2 > 9 8 5 ) IEO ..... .. . 

985 FOFMAT<5X**EON f *»I2) . . ... ....... 

VRITE<2>981) 

981 FOP PAT ( 2X* +TQT #*»2X»*PARAH »*»2X»*X VAR4»5X»4PRES$*»9X» 
SET UP DATA ARRAYS . 


♦TOT FP> 


IF(LATOP.EO.l) GO TO 800 . _. .. . ... 

DO 804 I*1»NPTS .... 

IF ( I EO-2 ) 601* 602* 803 

601 XtN*I>-Y(I)-2*HPG/CRH0PSAREA4VEim*VEim>*AXm 

GO TO 804 • 

802 X(N»I>«Y<I)"2*M*G/(RH0*SAREA4VEL(I)*VEL(1))*AZ(I) . 

C-0 TO 804 

803 X(N* I)*Y< II-2*IY/{RH0*SAREA4CBAR*YEL<I)*VEtt I))P(ODOT(U 

i-nz-ix)/iYpp(H4R(i)-ixz/iYt(R(u*Rm-pa)*p<im 

804 CONTINUE 

GO TO 805 . . 

800 00 6C6 I *1* NPTS 

IF ( I EO-2 ) 807*806* 809 

807 X(K*I>-Ym»2*M*G/CRH0*SAREA*VEim*VEltin*AYm 
GO TO 806 

808 XtN, n-Ym-2*IX/(RH0+SAREA*8SPAN*VEL(U*VEL(in*(P00T(H 
U IY-IZ) / 1X40 (I )PRm-{ IX Z/ IX I *<P(I >*0C I ) +RDOT { I > ) ) 

GO TP 606 

809 XtU,l>-Y(I)-2*IZ/<RHO*SAREA + 8SPAN*VEim*VEL ( 1 1 > ♦ ( ROOT ( I : 
u I x— I r > / 1 Z* P 1 1) *01 1>-( IXZ/IZ >*<PDOT( I)-0(I )*R< I) I ) 

806 CCNTINUE 

805 CONTINUE 

IANOVA-O $NM1*N-1 4N2M1«2*N-1 JIPASS-0 

N2-N2N14M2M1 SNAXL AG-NPTS/10 

LINCNT-0 SLMl-NHl 

F1-F2-FCRT 

IF ( I FL AG . EO. 1) LINOP-1 =.. 

DO 2C6 I * 1 * MAXI AG 
206 XL AG < I ) - 1— 1 

00 51 1*1*2401 
51 A ( I ) *0 • 

DO 52 I«l*Hf*l 


PROGRAM STEP 


FTN *.8*528 


82/10/28. 07. *1.22 


ACI»I+N)-1. _ 

00 53 1*1. NH1 

A(I+N>1)— 1. - 

DC 50 I*1»N 

XEAR(I)*XSUH{I)“ICNTtII*0. ... 

DO 50 II*1*N . . 

XXSUM! II»I)*0. 

00 ICO II*1.N _.. 

00 ICO I * 1 > N . 

00 100 J • 1. NPTS 

XXSUM! I, II )*XXSUHd» II)+X(X» 

on ioi II-i.n 

DO 101 3*1. NPTS .... . - - 

XSt!M(II)*XSUM(II)*X(II.J» 

DO 201 1*1. N 

DP 20C J * 1 » NPTS 

XPAR t I)*XPAP(I)*X<I»J) 

XBAP(I)*X8AP(n/NPTS 

DC 202 1*1. N 

DO 2C2 J-l.N 

SCI, J) -XXSUM (I. J>-XSUNm*XSUHCJ J/NPTS 
CONTINUE 
DO 203 I-l.N 
i SIGPAtI)*SCPT(StI.in 
DO 2C4 I-1.HM1 

1 PI - I *1 

DO ZCA J-IP1.N 

RRtl.Jl-SCI.Jt/JSIG-AdlTSIGMAU)) 

> RRt J,n*RR(I. J) 

DO 205 I-l.N 

s pr( x.n-i. ... ... ... . 

or 210 1 * 1 . N 

SIGNAC n-SIGHA«I>/SQRT(FLOATINPTS) S ... 
) CONTINUE 

PHI*NPTS-1 - ■ 

DO 3C1 1*1. N*U 

t 8(1 ) *0. 

SY* SIGMA (N)*SORT(RRCN»N) /PHI) 


REDEFINE XXSUM.X VITH REDUCED # OF DATA PTS FOR USE BY PRESS 


IF! IPRESOP.EQ.O) GO TO *50 




PROGRAM STEP 74/74 GPT*1~’.._ ’ * FTN 4*8 + 528 82/10/28. 07* 41.22 >AG£ ~ 4 


IF ( I PS KP .GT. 0) CALI PSET 
IF ( IPSKP .EO* 0) IPPTS-NPTS 

START LARGE L00° 

VM A X CALCULATED 


£ 145 


I-I * VMAX-O. SI PASS "I PASS + l __ 
IF(A(I#I )*GT*TOL) 250# 3C0 
IF ( ICNT ( I ) • EO *1 ) GO TO 300 
V<I)-A(I,N)*A(N#I)/ACI#n _ 

IF(V(I).GT.VMAX ) 260,300 
VM A X * V ( I ) SNHAX-I 

IFCLINOP.EO.l) LN1-LINMAX 

I F ( I • E 0* LM1 ) 330# 310 
I-1U 

GO TO 320 „ 

CONTINUE 
I-hMAX 

CALCULATE F 

IFtVMX.LT.TOL) 2000 #443 _ 

IF(LINCP.EO.O) 1 999# 2111 
F«PHI*VMAX/ 1 A(N# N)— VHAX) 

IF { F • GT *0# ) GO TO 444 
A(N#N)«VMAX SF*-F SNGO-1 SPRINT 998 
PRINT ^50# F# I ■ 

I F ( F *GT • FI . OR • L INOP* EO* 1 ) 400#1999 
IF ( I PR ES OP • EQ • 0 ) GO TO 403 

CALC PRESS 


CALL PRESS _ 

CONTINUE 

UPDATE THE A MATRIX __ 

ICMT(I)*1 SPHI-PHI-1 
DO 4C1 II«1#N2*1 

DO 401 JJ-1#N2M1 _ 

X F C 1 1 • NE • I ) GO TO 402 

AP(II# JJ)-A(II, JJ)/A(II#II) _ _ 

GO TO 401 

AP(II# JJ)»A(II# JJ)-A(II#I)*ACI#J4)/A(I#I) 

CONTINUE 

00 All 12 "1 # N2 Ml 


o > 

20 r*. 


O u 

-cr>- 

-F=-f>S- 





0X2 




CIIJiNOX+AI-M 6SV 
N'I-1I fcEV 00 
O-AI 

O-VAONVI 2EV 
3ivodn nvo 


s»oivo ssaaa ni oasn (a<»v> xxaivw v aivodn 


_ _ _ 2ev 01 os <0*03*d0S3adi>di 

3ANUN00 IEV 

' _ _ ' _ ' " OEV 01 00 

! ' _ ' - - - I» VAJNV1 

' ' ' lavwwns VAUNV M3N 

'■ ' ' 3r)NUU39 62V 

■ * ' ' <2r'2I)dV-(2f'2I> V 92V 

i IW2N'I*2f 9 2V QQ 

71 ! IH2N*I»2I 92*? 00 

30NI1N03 82V 

_ __ (Niux*N'Niwi*H)v/(rr^Ni»a)v*(NiHi*«‘ii)v-(f r‘n )v«< rr'inov i2v 
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S$DR-1.-A<N,N) SXNS0R-$$0R/IV ■ 

S$REG*A(N,N) SXNSSREG-SSREG/PHI 

PVAP-100+SSDR ■ ■ 

PRINT 954, PVAR - 

SCR FS*SGRT ( A < N> N) *S (N#N ) /PHI ) .... _ - 

PRINT 955, SORES • ••• 

FT0T-XNS0R/XMS3REG 

PRINT 957,FT0T ... - - - - 

FCPNAT<1X,*T0TAL F VALUE IS*»E12.3) 

PRINT SUHMAPY ON TAPE2 

VPITE<2»982) IV»NHAX,PVAR,PPRT(NNAXJ>FTOT 

F0RNAT(lX>15>2X,I5>3X>F6.2,3X»E12.5»3X»€12.3i .. 

TN( I V ) *FL OAT ( IV) : 

PPLT(IV)»PPRT(NNAX) 

FPLT(IV)-FTOT -■ 

OP 449 II-1.NN1 

IF( ICNTtlIJ.NE.il GO TO 448 ... 

B(II)*A(II>N) + SORT<$(N,NJ/S<II,in) ,~ 

STOEPni)*SORES*SQRTCA(Il»Il)/S(lI»ID) 

GO TO 449 - 

I B ( I I ) *0. *STOEP(II)»0. 

i CONTINUE 

Sl'N *0. 

DP 451 Il-1,NM1 ... . 

. SUH»SUN+B(II)*XBARCII1 

B(N )»XBAR(N)-SUM — - 

PPINT958 . 

PRINT 956» ( B( 1 1 )» II*1» N) __ 

PRINT 956, (ST0ER(I1)»II*1aNH1) 

DO 460 11*1, NPTS 

YHATCII)-B(N) — 

DO 461 JJ-1,NH1 , ... - .— 

L YHAT(in»YHAT(II> + B(4J)*X(JJ,Il) ....... 

} XNUai)*Y(II)-YHAT(II) . .. 

CAIL P AUTO t XNU> M AXL A gJnPTS^W# XXSUH(N»N)> IV* 1» YHAT# Y) " 
CALL INFOPLT<0»NPTS»T>l»Y»l,0.»0.>0.»0.»l.,4»4HTINt,t 
16HY,YHAT,22,5.,3.5,.75,.75) 

CALL INF0PLTC9,NPTS»T,1»YHAT»1,0.»0.»0.,0.#1.»4»4HTX 

16HY, YHAT, C, 5. ,3 .5, .75, .75) 



CALI INFOPLT<9>NPTS,T»1»XNU,1»0.»0.»0.,0.,0.,4»4HTIHE>6» 
16HT- THAT, 22, 5., 3. 5, 0.0, 4. 75) 

CALL INFOPLT<l,KAXLAG,XLAG,l,W,l,0.,0.,0.,0.,0.,3, 

13HLAG, 6, 6HAUT0 C» 22, 5. » 3. 5>0.C> 4.75 > 

. 470 CONTINUE 

IF(NGO.EO.l) GO TO 1000 

GC TO 2112 - - 

2111 LINOP-0 $L INCNT «L INNAX JLM1»N-1 ' 

GO TO 450 ' ' ' - 

2112 L INCNT *L INCNT + 1 

IFUINCNT.GE.LINNAX) LINOP-O 

LMi-N-1 ' 

GO TO 450 * 

950 FCRMAT(1X,///,10X, .MAXIMUM F VALUE IS *»E10»3»* FOR VARIABLE 

951 FCRMATUX,. VARIABLES IN REGRESSI0N*,24I4) U 

952 FCRMATdX, .PARTIAL F VALUE FOR VARIABLE *,I3,* IS *»E10.3) 

953 F0RHAT(1X».VARIABLE *,I3,* ELIMINATED*) 

954 FCRMATdX, ‘PERCENT VARIATION EXPLAINED IS *,F6.2) 

.. 955 F0PHAT(1X,.STD, DEVIATION OF RESIDUALS IS *,E10.3) 

. 958 FCRMATdX, .NEW PARAMETER ESTIMATES AND STD. DEV. ARE*) 

956 FORMAT dX,15E9.3,/,lX,10E9.3) 

998 FORMAT (IX. ‘NEGATIVE F VALUE CALCULATED*) 

1999 CONTINUE 

IF(IPRESnP.EO.C) GO TO 1000 ' ' 

C * L1 - p L0T2(TN,PPLT»FPlT» I V, 24HT0T4L NO. PARAM IN M0DEL#24, 

* 5HPPESS,5,4HFT0T,4) 

1000 CONTINUE - 

ENDFILE 2 
REWIND 2 
WRITE ( 6» 3C00 ) 

3000 F0PMAT(1H1,2X>* *) - 

DO 1500 1-1,100 ~ - - - - . 

RE AO ( 2, 2001 ) TITLE 

IF (EOF (2) ) 2500,1200 ' ' ' 

1200 CONTINUE ■ 

WRITE (6,2001) TITLE ~ ' 

1500 CONTINUE 

2001 FrP*'AT(8AlC) "" 

2500 CONTINUE _ 

GO TO 2502 ■“ ‘ 

2501 CALL CALPLT(0.,0.,999) 

STOP 
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SUBROUTINE DATAS ET ( TS,T» Y, X) ... - - 

REAL M,IY,IX,IZ,IXZ 

DIMENSION T(500 ) , Y < 500), X( 25, 500) 

DIMENSION BETA(500),ALPH(500) 

DIMENSION CEL A( 500) > DELR (500) ,DELE<500) 

DIMENSION NAMES ( 50 ) , IUNITS ( 50 )»HDR ( 8 ) ,DATAI50) 

DIMENSION PDO(500),ODP(500)»RDD(500)»PT£N(500)»QTEH(500),RTEnl5 

1*WRK1( 500) > WRY 2 ( 5C0 ) ,WRK3( 500) » VRK4( 500) »WRK5 ( 500) , WRK6 (500) 
COMMON/ AC DAT A/ S»B,C>M>P HO »G,IX,IY>IZ>IXZ, DELET, ALPHT 
1»PETT,OELAT»05LRT>OT>PT>RT 

COtfPON/FLAGS/ lPSKP»NPTS#IOIH»JDIN#NHAX#IHIN#ICNTt25)#I0RD(23) 

♦ , IPPTS>LATnPf ITRI^OP#TCALl> IACEIOP, IF HOP 

COHf*OM/OROEP/ IBQtH 

COMMON/ AC CEl/ AX ( 5 0 D > > A Y { 500 ) »AZ(500)# PD0T(50O)# QDQT(500)#RDQT( 
1 )*VEl< SCO)# P< 500 ># 0(500) #P (500) 


XF(ICAU.GT*l) GO TO 4»6 - 

J2-NPTS/2 
J PTS *2* J2 

I F ( CNPTS-*JPT$ ) .60* 0) NPT$*JPTS-X 

J DI **NPTS +N 

ID I H*N 

PE WIND 1 ... 

I D r 1 - 

NCH«20 

► READd) ID»NCH, C NAMES (I ), I“1»NCH) » ( IUNITS 1 1 )>I*1»NCH)»H0R 

t IF(E0F(1)> 999*, 502 

502 PEAOUJ (OATAU >, J*1,NCH) .... -- 

IF* EOF Cl) * 9996,8001 _ 

8001 IFC ITS). GT. DAT A(l) ) 502,600 
600 CONTINUE 

IFdTRIMCP.EO.l ) GO TO 602 

8ETT*DATA (3 ) t DE L AT* 0 AT A (13) *0ELRT-0ATA<15 ) JOELET -DATA ( 1 A ) 
AL P HT* DAT A( 19 ) 

602 PRINT 1980, ALPHT, B ETT, DE L AT, DELET, DELRT - 

I960 FCRMATdX,///,I0X,*TRIM VALUES*/, 15 X, 

♦♦ALPHT BETT AILT DELET DELRT*/, 

♦10X,5E12.5) 

DO 15 I*1,NPTS 

READ (I) ( DAT A { J ) , J*l, NCH) 
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~ * IF (EOF (1) ) 9996,601 

601 CONTINUE 

45 T ( I ) » OATAC 1 )-TS 

VE L ( I ) * DATA (2 ) _ " " 

BET A ( I ) » DATA ( 1 8 )-BETT ' 

ALPH(I)- DATAQ9J-ALPHT 

__ PCI)- DATA ( 5 ) _ 1 

50 ‘ ‘J* Q 1 1 ) ■ D AT A < 6 } ^ 

RID- D A T A ( 7 ) ‘ ‘ " ' "" 

AX { I ) • DAT A ( 10 ) “ ' 

AY ( I ) - DAT A ( 11 ) 

__ Aim* D A T A ( 12 ) “ 

55 DE L A ( I ) ■ DATA (13 )-DE l AT ^ 1 _ 

_ OELR(I)- DATA { 1 5 l-DELRT _ I J 

DEL E < I ) -P ATAC 14 )-DELET 

PCDT ( I ) -DATA ( 16 ) S0D0T(I)*DATA(20) $RDOT CD-DATA ( 17) 

15 CONTINUE 

‘60 “ IF(IACEL0P«EQ«0) GO TO 46 

___ CALI SECDER(3,3,T,P,PDO,PTEM,NPTS,PO,P3,.5,.5,VRKl,VRK2,VRK3, 

0__ _ 1WRK4,WPK5,VRK6) 

i_ II "" CALL SEC0ER(3,3,T,O#ODD,QTENfNPTS# PO# P3# «5# • 5» VRK1# WRK2# WRK3» 

1WRK4, WRK5, WRK6) 

65 "" CALL $ECDEP(3,3,T,R,RDD,RTEM,NPT$,P0,P3#.5,.5,VRK1,VRKZ»VRK3# 

: _ 1WPK4,WRK5, WPK6) 

1_1__ 21_" " 00 * 5 I -1 > N PTS 

PDOTCI >-DEPSP(T< I)>T#P»NPT5#PDD*PTEN) 

{_ ODOT ( I ) -DERSP(T( I >,T,Q,NPTS, CDD,QTEH) 

“ "70 RCOT(I)*DEPSP(T(I),T,R,NPTS,RDO,RTEfl) 

[I J, 45 CONTINUE 

; 46 CONTINUE 

‘ _ IF ( IF I LOP • EQ.O ) GO TO 47 _ ‘ 

„ “““ CALI F It ( T, AT , NPTS ) 

- 75 _ 47 CONTINUE 

_ _r_ c 

1_ I F ( L AT OP • NE • 1 ) GO TO 803 

DO 800 I * 1, NPTS _ _ _ 

X ( 1 , I ) -BETA ( I ) _ __ 

80 X ( 2 > I ) • P ( I ) * B/ ( 2 • *VEL < I ) > * _ _ 

X(3,n-PCI)48/(2.4VELCI)) 

X ( 4, I ) -DE L A C I ) 

X ( 5 , I ) -DEI R C I ) 

X(6,I)-ALPHCI)*xa,I) 









SUBROUTINE DATASET 74/74 OPT«l FTN 4.8+528 82/10/28. 07.41.22 PACE 3 


: 85 X(7iI)«AtPHm*)t(2»I) 

X(8,I)-XC3>n*ALPH<I) 

' . X (9* I J »DEl A * I ) PAIPM ( I ) '• 

_ _ _ _ XtlO»I)»CELR(I)*ALPM(I) 

■ Xdl.n-AtPHtl J4ALPH(I)*X(1,I) 

90 _■ _ X(12»I)*ALPH(I)*AIP4(I)*X(2»I) ;? 

• x{i3.n*AiPHcn*AtPH(i>*x{3»i ) 

_ _ X(14*I)«ALPHm*ALPHm*X(4,I) 

__ X(15*I)«AlPH(IJ*AlPHtI)*X(5»I) • 

x(i6*n«xa.n*x(i»n J ~ 

.95 X(17*I)«X(1»IJ*X(16»I) ' 1" _ 

x(i8,i)»xa,n*xii7/n ' _ "" 1C1 

x{i9,i)-x(i,n*x(i8»i) _ 1 

; X(20*t)»X(l*I)*X(ll»I) * ■ ■' 

X(21,I)»ALrH<n*X(16,I> 

100 X ( 22* I ) »Al PH( I ) 

_ _ _ X(23.I)*ALPH(I)4ALPH(I) 

’ " X(24*I)»ALPH(I)*X(23*1) ^ 

w ~ 800 CONTINUE 

' - GO TO 804 _ ■ ~ 

" 105 _" 803 CONTINUE 

DO 8C1 I«l* NPTS " 

_ ' 1" X (1 * I ) »AL PH ( I ) 

_ _ “■ X(2,I)-C/(2*VFLCin*om 

; X ( 3 * I ) "DELE ( I ) J ~ ~ ~ ~ 

_ _ lio x(4,i).xa,n*x(i,i) 

__ x(5»i)«x(i,n*x(2»i) 

; _ X(6» I ) «X ( 1 » I ) * X ( 3* I ) 1 

X(7.I)*X(4,n*X(2#I) • _ 

_ _■ _ X(8,I)«Xt4,I)*Xt3*II *" 

115 X(9*I)«X(4*I)*X(1*I) 

X(10#I )*X(9,n4X(l,I) 

■ _ x(ii,ii-xtio>n*x(i»n 

: x(i2*i)«x(ii*i)*x(i»l) 

x(i3.n-x(i2,i)4x(i.n _ 1 

120 801 Xtl4»I>«X<13»II*X(l,n 

804 CONTINUE ; ' _ 

_ PRINT 964, ID. IEQ*NPTS 

IF(ICAU.GT.l) GO TO 999 
IF< LATOP.EO.O) GO TO 50 
125 PPJNT 968 

968 FCPNAT ( 1 X * 7 X , »TIHE*,UX,*V*,12X,*BETA4*11X,«P*,14X,*R*»11X* 


j 
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SUBROUTINE DATASET 
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PAGE 


. 13 ° 


135 


... 1*0 


^ 145" 


' 150 


1*DELA«.11X,*DSLR*> 

964 FCRMAT(1H1»3X»*RUN*»I3>3X»*EQUATI0N *» 12/ 3 X> *NPTS*> I3» /// ) 

PRINT 962»(T<I)»VEL(I)»BETA(I)»P(I)»R(I)»0ELA(I)jOELR(I)»I*l»NPTS) 
<562 F0P.NAT(7(2X,E12.4) > 

CALL INF0PLT(1>NPTS»T>1>DELA*1>0.»0.»0.»0.>0.»1»1HTj 
14»4HDELA»22»7.»3.».73».75) 

CALL INFCPLT<l»NPTS>T»l.DELR»l»0.>0.fO./0.»0.»l»lHT» 
14,4HDELR»22.7.» 5.. .75, .75) 

GO TO <599 

30 CONTINUE 

PRINT 998 

998 F0RNAT<1X,7X>*TINE*#11X>*V*>11X**ALPHA**11X>*0*>14X#*0ELE*) . 

PRINT 997, <T<I)> VEL(I), ALPH( I) » 0 (I) , DELE ( I ) » I*1»NPTS> 

997 FOPNAT ( 5(2X>E12.4) ) . 

CALL INFnPLT(l>NPT$>T>l>DELE»l>0»>C»>0»,0.,0«»l>lHT» . 

14,4H0ELE,22,7.,5.,.75».75J 

C-0 Tn 999 

9996 PRINT 9999 

9999 FORNAT (1X,*E0F ON DUNNY READ*) . ... 

GO TG 999 

9996 PRINT 9997* J,I . 

9997 F CP NAT ( IX » *EOF ON DATA REAO*, 5X» *INDEX J>I» *,2110) _ 

C-0 TO 999 

9994 PRINT 9995 

9995 FOP NAT < 1 X, *EOF ON HOR PEAO*) ... . 

999 CONTINUE . 

RETURN ....... 

END . 
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7*/74 OPT«l FTN 4. 8 + 528 

FUNCTION DERSP<XX,X,Y,N,P,H) _ _ 

FUNCTION DERSP ' " “ — ' 

DERSP IS USED TO OBTAIN THE FIRST DERIVATIVE OF SPLINE 
CURVE FITTED DATA 

USAGE - 

X • DERSP(XX,X»Y,N,P,H) ' 

NOTE - IF XX LESS THAN X(l> THEN DERSP - DX(1)/DY 

IF XX GREATER THAN X(N) THEN DERSP - DXCNI/DY ' 

VHERE - -- ' 

xx INDEPENDENT VARIABLE FOR WHICH INTERPOLATED SLOPE 

IS DERSIRED 

X N-DIMENSIONEO VECTOR OF INDEPENDENT POINTS 

Y N-OIMENSIONED VECTOR OF DEPENDENT POINTS 

N NUH9SR OF DATA POINTS 

P N-DINENSIONEO VECTOR FROM UPDATE 

. H (N-U-DIHENSIONED VECTOR FROM UPDATE 

SUBROUTINES CALLED - ' • 

NONE 

OIHENS ION X(l),Y(l)»Ptl),Hm 

XP-XX 

IFIXX.LT. Xlll) GO TO 1 

K-N-l 

DO 2 I-1,K 

IFIXX.LT. xa+in GO TO 3 “■ ■ ■ 

CONTINUE 

I *K ~ — ’ 

XP-X(N) J 

GO TO 3 ' ' * ‘ 

xp*x(D ' : ~ - - — 

F1*(X(I+IJ— XP)**2 - 

F2* ( X P-X ( I ) )t*2 * 

F3-HCI )/3. ‘ 

rERSP-((F3-Fi/H(in*pm ♦ <F2/H<n~F3)*pn+i>)/2.+tY(i+i>-Ym)/ 

L H ( I ) 

RETURN 

ENO 


82/10/28. 

DSCF 

230 

DSC F 

10 

DSCF 

20 

DSCF 

30 

DSCF 

40 

DSCF 

50 

DSCF 

60 

DSCF 

70 

DSCF 

eo 

DSCF 

90 

DSCF 

100 

DSCF 

110 

DSCF 

120 

DSCF 

130 

DSCF 

1 AO 

DSCF 

150 

DSCF 

160 

DSCF 

170 

DSCF 

180 

DSCF 

190 

DSCF 

2 CO 

DSCF 

210 

DSCF 

220 

DSCF 

2 AO 

DSCF 

250 

DSCF 

260 

DSCF 

270 

DSCF 

280 

DSCF 

290 

DSCF 

300 

DSCF 

310 

DSCF 

320 

DSCF 

330 

DSCF 

3 AO 

DSCF 

350 

DSCF 

360 

DSCF 

370 

DSCF 

380 

DSCF 

390 

DSCF 

400 

DSCF 

410 

DSCF 

420 
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SUBROUTINE $eCDFR(llH2>X,Y*P*Hj>N,P0,P3,XKl,XK2,A,B,CfD,GAHNA, 
1 BETA) 

SUBROUTINE SEC DER 

SEC DE R IS US E 0 WITH FUNCTION SPLINE TO PREFORM A SPLINE 
INTERPOLATION. IT IS USED TO GENERATE P AND H. 


USAGE - 

CALL SECDER(LljL2*X>Y>P#H#N*P0#P3>XKl>XK2#A#B>C#0>GAMMAfBETA) 
WHERE “ 

Ll» L2 OETE RHINE THE END CONDITIONS AT X (1 > AND X(N) TO BE 
USED. (SEE BELOW) 

X N-DI MENS IONED VECTOR OF INDEPENT POINTS 

Y N-P I *ENS I ONE 0 VECTOR OF DEPENDENT POINTS 

P N-OIMENSIONED VECTOR TO BE RETURNED 

H (N-l )-DIHEN$IQNED VECTOR TO BE RETURNED 

N NUMBER OF DATA POINTS 

SECOND DERIVATIVES ARE GIVEN AT THE END POINTS 
XK1 NOT US EO 

XK2 NOT USED 

IF Ll-1 THEN 

PO SECOND DERIVATIVE AT XU). Yd) 

IF L2-1 THEN 

P3 SECOND DERIVATIVE AT X(N),Y(N) 

FIRST DERIVATIVES ARE GIVEN AT THE END POINT 
XK1 NOT USED 

XK2 NOT USED 

IF U-2 THEN 

PO FIRST DERIVATIVE AT X(1)»Y(1) 

IF 12*2 THEN 

P3 FIRST DERIVATIVE AT X(N),Y(N> 

NO INF OPMAT I ON ABOUT THE CURVE IS KNOWN 

PO NOT USED 

P3 NOT USED 

IF Ll-3 THEN 

XK1 P'*(3»0) - XK1*P* • (3.1) » XK1 GREATER THAN 0 . 

IF L2*3 THEN 

XK2 P • * ( 3 » N ) « XK 2* P' ' ( 3. N-l ) » XK2 GREATER THAN 0 
A N-OIMENSIONED WORK VECTOR 

B N-DI MENSIQNFD WORK VECTOR 

C N-D I MENS I ONE D WORK VECTOR 


UPO 

480 

UPO 

490 

UPO 

10 

UPO 

20 

UPD 

30 

UPD 

40 

UPD 

50 

UPD 

60 

UPO 

70 

UPD 

80 

UPD 

90 

UPD 

100 

UPD 

110 

UPD 

120 

UPO 

130 

UPD 

140 

UPD 

150 

UPD 

160 

UPD 

170 

UPD 

ieo 

UPD 

190 

UPD 

200 

UPO 

210 

UPD 

2 20 

UPD 

230 

UPD 

240 

UPD 

250 

UPD 

260 

UPD 

270 

UPO 

280 

UPD 

290 

UPD 

300 

UPD 

310 

UPD 

320 

UPD 

330 

UPD 

340 

UPD 

350 

UPD 

360 

UPD 

370 

UPD 

380 

UPD 

390 

UPD 

4C0 


- © o 

-TjS 

o g- 
-o-S- 

20 jp* 

' <0~ tf 
' a >- 

> Q 

r° m - 





0 N-DIHENSIONED WORK VECTOR 

BETA N-OIMENSIONEO WORK VECTOR ._ 

GAMMA N-D I HEN STONED WORK VECTOR 

SUBROUTINES CALI - .......... ... 

NONE 

DIMENSION X(N)»Y(N)» A l N ) »8 (N ) » C ( N ) , D (N ) »GAMH A IN 
K-N-l 

00 1 J-1»K 

1 H ( J ) *X C J+D-X ( J ) 

DC 2 J *2> K 

A ( J ) • . . ...... ... . 

B(J) • 2.*<H<J)+H(J-1)>/HIJ> 

C(J) * I. 

2 0 ( J ) - 6./H(J)R((Y<J+X)-Y(J)J/H(J)-(Y(J)-T<J-l) 
IFtLl.E0.2I GO TO 20 

IFCL1.E0.3) GO TO 10 

B<1>«1. .. . ... _. ._. 

C < 1 1 -0 . 

om-po 

GO TO 30 ... 

10 8 (1 )-l. - _. 

cm— xki .. 

D ( 1 ) *0 • .... 

GO TO 30 _ . . - 

20 8(1 

CIll'HIll/t. ... 

Dll>-(Y<2J-YUn/Hm-P0 

30 IF(L2.E0.2) GO TO 21 

IFIL2.E0.3) GOTO 11 

A(N )*0. .... 

B(N)«1. 

D(N )-P3 .. . 

GO TO AO . 

11 Atm — XK2 _. ... 

B<N)»1. . ... : 

D(N )*0« .. 

GO TO 40 . - . 

21 A(N)«H(K)/6. 

B(N)-H(K)/3. 

D(N).P3-(Y(N)-Y«KJ J/HtK) 


),BETA(N),H<N),P<N) 







POOR 




056 Gdfl 

o v 6 adn 

0£6 GdD 
026 QdO 
0X6 adfi 
006 adn 

068 QdH 
088 adn 
oi 9 adn 
099 QdH 
058 QdO 


. ' GN3 

Ndni3d 

(W)Vi3a/lt+W)d*<W)3-<W)VW^V0- lW)d l 

T-a-H 

x*i«r z. oa 

(N)VW^VO - JM)d 

(r>vi39/<<i-nvwwvo*cr)v-(na)-(r ivriwvo 9 

a-D U38/C W)D*(f >V-< r>9«( DVX39 

N*2»r 9 aa 

(I)Vi38/{I)a-(I>VW^IV0 

" ~ (t)9*(I)V138 OV 
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SUBROUTINE Fit 


74/74 CRT" 1 


FTN 4*8+928 


1 


5 


10 


15 


20 


25 


30 


35 


40 


_ SUBROUTINE FIL<T,P,NPT$> 

DIMENSION PF<500)> Ht 500) # P ( 500 ) # T ( 500) 

PI-3.14159 SFC-2.5 SFT-2.51 SWC-2*PI*FC SVT*2*PI+FT 

NNI0-NPTS/2 SOT". 05 
V2* (VT-WC ) WT-VC ) 

DO 3 I"1#NMID 
— ' K * X 

‘ 3 H(I ) *P 1/ (2*K*DT ) *(SIN( VT*K*OT )*SIN (WC*K*DT) ) /(PI+PI— V2*K*DT 
1*K* PT ) 

HO-FC+FT 
NPM1-NPTS-1 
HNDPH* HO 
DO 1 I "1# NNID 

1 HNORH" HNORK+H ( I ) *2 • 

DO 2 I "1# NNID • 

2 H < I ) *H C I ) /HNORM 

HC-HO/HNORN 
DO 5 I "2 » NN ID 
I HI * I— 1 

PFm-HC + P{I) " ' “ 

DO 51 J ■ 1 # I FI 1 “ * ” 

51 PF( I ) • PF ( I ) *H C J ) ♦ ( P( I+J } 4 P( I— J > > 

DO 52 J*I# NMIO 

52 PF(I)-PFCI)+2*H(J)+P(I+J) 

5 CONTINUE 

NP2 * NM ID + 2 ~ " 

DO 4 I«NP2#NPN1 _ 

NPHI-NPTS-I ' 

PF(I)»HO*P(I) * ~ 

DO 41 J*1*KPMI ' “ * 

_ 41 PF(I)-PF(I.)+H(J>«(-PtZ-J) + P<Z+Jir 

NPHIP1 *NPMI+1 
DO 42 J-NPPI PI# NMID 
42 PF( I > — PFfll *2*H t J ) +P( I— J ) 

4 CONTINUE 

PF (1 )-HG*.PCl) SPF(NPT$)*HO*P(NPT$) SPF ( NMID+l ) «HO*P (NHI DU ) 
DO 1C J-UNMID 
PF(1)"PF(1)+2*H{ J)*-PU + J) 

PFtNNIOU)*PF<NNIDU)+H< J > ♦ C P C NH IDU ♦ J > *P t NH I DU-4M 
10 PF(NPT$)«PF(NPTS)*2*H(J)*P{NPTS~J) 

DO 6 I"1#NPTS 

6 P<I>-PFm 





SUBROUTINE PLQT2 


PTM 4.84*528 


82 / 10 / 28 * 07 * 41*22 


PAGE 
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^ SUBROUTINE PL0T2(T>X>Y*N>lABT*LTjLABX>lX*tA8Y*lY) 

' DIMENSION T<1),X<1>*YU) 

Nl-N+l - . 

' " N2-N+2 . . 

5 CALL ASCALE(T>4.fN,lflO.) 

CALL A$CAtE(X,2.f N,lf 10.) 

- • • - CALI A$CALE(Y>2.,N,1,10*) . 

— ' CALL CALPLTd.fi. #-3) 

~ CAU AXES f0*>0.# O.f 4.fT(Nl)f T(N2)f l*f 5.f L ABT > • 13 f *-LT* 2) 

XO CALL CALPLT(C. .l.f-3) ... -- 

_ " ' CALI AXES(0.f0.#90.#2.#Y(Ni)#X(N2)Vl.f5.*LA8X#.15#LX#2) 

... CAlL xiNPLT(TfXfNflfOfOfl) . 

- ... - - CALl CALPLTtO.f 2.5>-3) : 

CALL AXES(C.fC.f90.f2.fY(Nl)fY(N2)fl.f3.fLABYf *15i*lYf2) 


X5 CALL LINPLT(TfY#Nf lfOfOf 1) . - . - - - - - -f 

CALL N FRAME _ . . _ J 

RETURN . - - i 
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COT 

( ( T ' I ) > ) 1 ^ < s ( t )vn"TS 
fj O t a T r<i? nr! 

f i/d 


in t‘T i v u*3 

ci on / ( p i no cy * { i ) nncx -{ r* I) no <;v x= ( r M »9 

n* T = r <>oz .in 

T b r r*. 7 ni| 
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, ~i hn / ( n dx = ( n *v . 1 * 

T y ? 
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'J 0 J 
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■no m no 

wminiio 
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{ T ) 7 V *< ( T >11A*< n T3A* Vl>VS*( H>0 / C )*W*?«( M AWT^IV c « >* 

V- ui no 

c n vv *m nr’A*u ) i 3 a*v:« fvs*r*MH) /o*w*? =f i ) a» ( i 'n > x i* - 

0 id n' T * r b^vy nr* 
mO (.1 1 fio f T V) ‘ • dfi 1*0 ) -IT 


" A 7 J > ,w v J V r t jfi i >; 


f X 'A ' i 'S J ) 1 ' S V IV .1 11 VO 
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07 


co 

\D 


07 


f'V 


00 


'■'O 


0*7 
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F 1 
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I'J 


O' 

K) 


3 v/ 


3 b 


A J 


P- i;. - . jl'r L< 1 r.PjT, i; :-| P1J r , T m'M : 1 * I * MP ) f , T A P f 0 * l • U I FM.I I ) 

K 1 m L .,1Y, iX,ii,iXt 

P i .;.J I ! N Pk ( / 7 4 ) > A { 7 -it ?■-) ) , AP { M, 7 7 ) 

F ■ I >V N J I ■ r : r l yUD) , Y { . I J # Y I A T ( V ' J o ) # X : 1 IJ l '#uC ) 

I’] r.i IO ! ‘.N 3 ( to a v ) > M p) W l** 3 ) 

lJ I M i h > 1 ■■} « 3 I Ui> h { A C ) , t- P A tv T ( A 3 ) 

t ■ i H M li, f) M < )>:■ ) , X L Ab ( f v, ) 

PI h, r-o II‘N li QN (3 ),KNT( 23) 

L.,r.r.:n /\ I a vt / x( Au, vuj ), xx sum ( ) , xl>uM <♦. > , xuak mo ) » b Mr a ml) 

l, ir.-ili J // r f > AT A/ 5AJ- EA, l 3 .\k, ■-)> r.HU> M : X , 1 Y, i 2 i tXZ,UF.LH ,ALPnT 

l T I > u l T > Ul L * T , C f > P T > P T 

..Ml'-Ib i /i Li- GS / NP T S , I Cn F ( A 1 ) , L A Ti ii> , 17 R I N 0 P , I L AL L > I A Ci Lit P, IP I L HP 

/ i 'kii >’ / I : Q , N 

l ;ii j m /u tl lL / ax ( tfC-./ ) , ay i j ■.)■.' ) M / ( m. ) > pp m ( / j.; ) , oi j[ ( )#> luji i o 

1 )f VI L( -K b) ,P MCl ) , U14 J„) M IV )b ) 

CPNMUN /« N J 1 / XN LT ( 2 3 ) 

LOl' iVA Li fi«. 1'. (A,FP) 

NAN » L t :>T / IN PU T / T S ^ N . y > 1 . OM, 

« NPTb» IPl JT, 1 H L A '3 > 

^ 3 A t* t A t H i> P A N p C :i A K > 

* 

* I X, iY, i 1XZ, 

* DlL r T ,ALPHF, ^ - TT, O' IAT> D. T, JT# P T, RT, 

* hC-< TMTP i NUP> la TJP, 1ACK LbP? IP 1 lUP 

A l p b T * : i ( T T * ( • l L L T * I. L L f T * [) t l T * Q l * P T ■ F T * 0 * 

CALL P^lUi-Ij 

2f.<? Tf L = X. 0[ UCALL=v 

u { ' AQ li=i, L3 

a -J Xf< Tl I ( J i ) c # OM 72 b5 + t 1 i-i ) * * J 17 A ^ i 
kl AM!, , INPUT ) 

It (: Mt M ) ) 2 V;) 1,2503 
2 5 )J K- IT Lib, INPUT) 

OF! j U t=l, 23 

3 0 k.n T (1) = lriX(XNL)rtl ) * b 7 • 3 ) 

PP Pit 45 4, (1,1*1, 17), iKNfi l),l = i»l7) 

N^c , l \ MA X- J 

i f ( l a r i tp * b y • i ) ).* 2 3 

IP ( L A! i)p • y. 1) L I N PA X = 5 

i>0 10 j J L*i,flbO 

i . l C * It f JN ( L ) Tf: Uu 3 v 

K A L L = 1C AL L+ 1 



APPENDIX 5 


This appendix contains the listing for PROGRAM STEPSPL (as used in 
example 2), 
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A2VIPJP. 82/10/28.NASA/LPC CY175-R NOS 1*4 <R 23 j 


0 7* 

0 7* 

C 7. 

07 

07. 

C7 • 

07. 

07* 

07. 

07. 

0 7. 

07. 

C7. 

07. 

0 7. 

C 7 • 
07. 
07. 
07. 
07. 
O' 0 7. 
° 07. 
07. 
07. 
07. 
07. 
07. 
07. 
07. 
07 
07 
07, 

0 7, 
07. 
07. 
07. 

C 7. 
07. 
07. 
07. 
07. 
07. 
07. 
07. 
07. 

07 
i 07 


41*20. STEPJ# T3000. 

*1,20. BATTER SON 

*1.20. USER, 043A50N. ' - 

*1. 21. CHARGE, 101216,LRC. 

*1.21. GET, STEP1. 

4 c 2 . F TN ( I » S TE PI , R -0, A# P L *1 5000 ) 
*1.*32.‘gET,0UH. 535 CP SECCNDS COMPILATION TIKE 

41.33. 'i£T#BlHVDP. • -• 

4 1.3 4. 5K IP?, SINVDP#7# # B. 

41.34. C0PYBR* PINVDP,TAPE1,1 # 

41.34. COPY COMPLETE. 

41. 34. REWIND# TAPE1. 

41. 34. ATTACH, FTNP-LIB/UN-LIBRARY. 

41. 34. 4TTACH#AKCLIB/ UN -LIBRARY. — 

4 1.3 5. ATTACH, L PC GO SF /UN- LI 3R APT. " 

41. 35. GET, ISSILIR/UN-474750C< 

L. 1 1 A. _ 


NOS 1.4 531 R CPFS 


R2142 B 


*; I • 35. G E T# I S S I L I R/ UN-474750C « - - ..... 

M:?6:LGcfr T >5K! e ' FTN " LI8/AKCLIB/LRCG ' OS,: ' ISSILI8 » PR ESETA.NGINF,KAP.SBEX) 

41.42. STOP 

41.42. 251100 MAXIMUM EXECUTION FL. ‘ 

^^•42- 8.700 CP SECONDS EXECUTION TIME. 

41.42. P10T*VARIAN<FPAMCNT»30#EDIT(1#30)) 

4-1.44. VO 01 

41.53. 10 FF A M ES / 2.31 METEPS GFNFRATFD 

*l*^ ,P J*I URE IKAGE FUE WIU BE saved ON OISK 

l 1 -:*'** PLCT 0UTPUT COMPLETED ..... 

^ , u^I* VARIAN ' (FRAf,CNT ' 30 ' EDIT(31 ' 6 0>) ~ ' - 

L.5b.VC02 • ■ • ■ .... ..... 

A2.C4. NO PLOTTING ATTEMPTED - 

M*^*w!: CT * VARIAN(FRAMCNT * 3 0'EDIT(61,90n 

42.C7.VC02 

42.13. NO PLOTTING ATTEMPTED 

4L # l^VC02 #VARIAN(FRAWCNT ’ 30,EDIT(91 ' 120n 
42.21. NO PLOTTING ATTEMPTED 

42^24^002 * VARlAN(FRArCHT " 30,EDITa21,:I50n 
NO PLOTTING ATTEMPTED 


42.30. 

42. 30. EXIT. 

42.30. UEAD# 

42.30. U EPF# 
42.70.UEMS# 


0 *002 KUN S • 

1 .426KUNS • 

21. 4 62K UN S . 

42.30. 175 CPU SEC - UECP/5.0* 

42.30. UECP# 43.379SECS. 

42.30. AESR, 192.282HNTS. 

42.30. APPP 0 IJM ATE JOB EXECUTION COST 



APPENDIX A 


This appendix contains a sample procedure file (JCL deck) for running 
example 1 at the Langley Research Center computer complex. 
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STEPVISE REGRESSION SUMMARY 


EQN * 1 

_JOT # P A R Art # X VAR 
”"l 1 93.43 

2 3 100.00 

EON $ Z 

TOT C PARArt f X VAR 
1 2 58.76 

1 __ 2 _ 3 72 . 6 * 

3 1 100.00 

EON f 3 

TOT # PARArt # X VAR 

1 2 BO • 54 

2 3 97.22 

" 3 1 99.94 

_ 4 ? 99.97 


PRESS 

-R 

TOT F 

' • 85710E + 03 


-R 

. 59773E+10 

PRESS 

-R 

TOT F 

• 53427E + 02 


-R 

• 53111 E+OZ 


— R 

*2 6 453E + 11 

PRESS 

•R 

TOT F 

• 1 6970 E 403 


-R 

.69937E+G3 


-R 

• 20 8 66 E + 05 


-R 

.27162E+05 


Ln 




PERCENT VARIATION EXPLAINED IS 99. 94 
STD. DEVIATION OF RESIDUALS IS .X48E-02 
TOTAL F VALUE IS .20S66E+05 
NEW PARAMETER ESTIMATES AND STD. DEV. ARE 

-.105E+01 . 153E+02-.996E+000. 0. 0. O.O.O.O. 0 , 0. 0. 0* -.731E-01 

' .2366-01 .8fc6E-CX .7846-020. 0. 0. 0. 0» " 0. _ 0. 0. 0. 0. 0. 

RESIDUAL MEAN IS -.52800E-14 
SIGMA SO OF RESIDUAL IS .832856-04 


MAXIMUM F VALUE IS .305E+02 FOR VARIABLE 7 _ _ 

VARIABLES IN REGRESSION lllOOOXOO'O OOOO _ __ __ 

PARTIAL F VALUE FOR VARIABLE l IS .296E+04 ~ " ' ----- • ^ J 

PARTIAL F VALUE FOR VARIABLE 2 IS .333E+04 _ _ _■ ' _ 

PARTIAL F VALUE FOP VARIABLE 3 IS .2O3E+05 

PARTIAL F VALUE FOR VARIABLE 7 IS .297E+02 

VASIA3LES IN REGRESSION 111C0010000000 
PARTIAL F VALUE FOR VARIABLE 1 IS .296E+04 

PARTIAL F VALUE FOP VARIABLE 2 IS .553E+04 

PARTIAL F VALUE FOR VARIABLE 3 IS .203E+05 __ _ _ _ _ 

P A P T I AL F VALUE FOR VARIABLE 7 IS .297E + 02 ’ _ 

PERCENT VARIATION EXPLAINED IS 99.97 

STD. DEVIATION OF RESIDUALS IS .1X2E-02 _ _ • 

TOTAL F VALUE IS .27162E+05 

NEW PARAMETER ESTIMATES AND STO. DEV. ARE 

-.1G7E+01 .143E»02-.991E+CC0. 0. 0. «327E*030. 0. 0. 0. 0... 0. ...... 0. ... -•7176-01 


•179E-01 .6 57 £-01 .595E-02C. 0. 0. 

RESITUAL MEAN IS -.55C55E-X* 

SIGMA SO OF RESIDUAL IS .47891E-04 


• 203E+020. 0. 0. .. 0. . 0. . 0. 






MAXIMUM F VALUE IS 


424E-A01 FOR VARIABLE 14 



EQUATION 3 


NPTS 43 


RUN 1 


WAX I NUN F VALUE IS .17AE+03 FOR VARIABLE 2 
VARIABLES IN REGRESSION 0 1 O 0 0 0 00000 

PERCENT VARIATION EXPLAINED IS 80. 

STD. DEVIATION OF RESIDUALS IS .25SE-01 ~ 

TOTAL F VALUE IS .16970E+03 

NEW PARAMETER ESTIMATES AND STO. DEV. ARE ' 

0* .19SE+020. 0. 0. 0. 0. 0» 

0. .IA9E + C10. 0. 0. 0» 0. 0* 

RESIDUAL MEAN IS -.A9282E-1A 

SIGMA SO OF RESIDUAL IS .2665AE-01 


Ln. MAXIMUM F VALUE IS .2A6E+03 FOR VARIABLE 3 

o> VARI A9L6S IN REGRESSION 0 1 X 0 0 0 0*0 ' ~0 0 " 0 

PARTIAL F VALUE FOR VARIABLE 2 IS .686E+03 

PARTIAL F VALUE FOR VAPIABLE 3 IS .240E+03 

VARIABLES IN REGRESSION 0 1 1 0 0 0 0 0000 

PARTIAL F VALUE FOR VARIABLE 2 IS .6S6E+03 

PARTIAL F VALUE for VARIABLE 3 IS .240E+03 

PERCENT VARIATION EXPLAINED IS 97.22 

STD. DEVIATION OF RESIDUALS IS .976E-02 

TOTAL E VALUE IS .69957E+03 

NEW PARAMETER ESTIMATES AND STD. DEV. ARE 

0. .160E+C2-.85SE+CCC. 0. 0. 0. 0. 

C. .S71E+C0 .517E-C1G. 0. 0. 0. 0. “ 

RESIDUAL MEAN IS -.73430E-14 

SIGMA SO OF RESIDUAL IS .38071E-02 


MAXIMUM F VALUE IS .173E+0* FOR VARIABLE ~ X ~ 

VARIABLES IN REGRESSION 11X00000000 
PARTIAL F VALUE FOR VARIABLE 1 IS .170E+0A 

PARTIAL F VALUE FOR VARIABLE 2 IS .260E+05 

PARTIAL F VALUE FOR VAPIABLE 3 IS .I21E+05 

VARIABLES IN REGRESSION 11100000000 
PARTIAL F VALUE FOR VARIABLE 1 IS .170E+04 

PARTIAL F VALUE FOR VARIABLE 2 IS .260E+03 

PARTIAL F VALUE FOR VAPIABLE 3 IS . 121S+05 



0 0 




:0+3T2T*- 



ot-avci‘5** si i?noiS3a du os vwois 
et-3a* 96A*- si Nvaw ivnoisaa 

•0 •OiO-3509* *0-3899* VO-3281* 


0 



•0 * 0T0+3CII *-20+ 3003* -10+3 00£* _ 

3aV • A30 * QiS ONV SBlVriliSa «313«Va7d «3N 
II+3ESV82* SI 3m»A d 1V1D1 
£0-3VIt* SI ST?nQIS3a dO NQI1TIA30 *ois 
00*001 SI Q3NI»ldX3 NOIiVla»A lN3D33d 




O OC 


RUN 


EQUATION 2 


NPTS 43 


MAXI FU* F VALUE IS •599E+02 FOR VARIABLE 
VARIABLES IN REGRESSION 0 1 0 0000 

PERCENT VARIATION EXPLAINED IS 58*76 
RTn. DEVIATION OF RESIDUALS IS .334E-01 
TOTAL F VALUE IS *58427E+02 

NEV PARAMETER ESTIMATES AND STD* DEV* ARE . 

0# -• 15CE +020 • 0. 0. 0« --- - 


. * 1 9 6E *0 10 • 0* C 

ESI DUAL ME AN I S -*79977E-13 
SIGMA SQ OF RESIDUAL IS .45829E-01 


_ 0* 0* _ 


2 

0 0 


0 *. 

0* 


0* 

0* 


0* 

0. 


— •131E+01 


Ln 


VAR I 
PART 
PART 
V ARI 
PART 
PI PT 
P E R C 
STD. 
TOTA 
NEW 
0* 


MAXIMUM F VALUE IS 
A3LES IN REGRESSION 0 
IAl F VALUE FOP VARIABLE 

IAL F VALUE FOP VARIABLE 

ABIES IN REGRESSION 0 
IAL F VALUE FOP VARIABLE 

IAL F VALUE FOR VARIABLE 

FNI VAP.IATICN EXPLAINED IS 72.64 
DEVIATION OF RESIDUALS IS *276E-01 
L F VALUE IS • 53 11 1 E +02 
PARAMETER ESTIMATES AND STD* DEV* APE 
-* 178E ♦C2-*705E+OC0 • 0* 


•208E+02 FOR VARIABLE 
l 1 0 0 o 0 

• 105E + 03 

• 203E+02 
0 0 0 0 

• 105E +0 3 

, 2 03E + 02 


IS 

IS 

1 

IS 

IS 


04 *161 E+01 .146E+CC0. 0* 

RESIDUAL MEAN IS 8 5 596E-1 3 

SIGMA SO OF RESIDUAL IS .30403E-01 


3 

0 0 


0 o 

0 0 


0* 


0* 

0. 


0. 

0* 


0 0 
0 0 


0 

0 


0* 

0. 


0. 

0* 


0* 

0« 


0* 

0. 


0* 

0< 


. . -i 

1 

-.13VE+01 


maximum f value is 

VARIABLES IN REGRESSION 1 
PARTIAL F VALLE FOR VARIABLE 
VALUE FOR VARIABLE 
VALUE FOR VARIABLE 
IN REGRESSION 1 
VALUE FOR VARIABLE 
VALUE FOP VARIABLE 
VALUE F OR_V ARI ABLE. 


PARTIAL ? 

partial f 

VARIABLES 

"partial f 
partial f 
‘partial f 


.239E+11 

for variable 

1 

1 

0 0 0 

0 

1 

IS 

•234E+11 


2 

IS 

.749E+11 


3 

IS 

.299E+11 


1 

1 

0 0 0 

0 

1 

IS 

.234E+11 


2 

IS 

.7V9E+11 


3 

IS 

. 2V9E + 11 



1 

0 



MAXIMUM F VALUE IS .B78E+03 FOR VARIABLE 1 
VARIABLES IN REGRESSION lOOOOOOOOOOOOO 
PERCENT VARIATION EXPLAINED IS 95. A3 

STD. DEVIATION OF RESIDUALS IS .139E-02 -■ 

TOTAL F VALUE IS .85710E+03 " • " ' ’ ' 

NEW PARAMETER ESTIMATES AND STD. DEV. ARE / 

.65 1E+0C0. 0. 0. 0. 0. . . 0 . 0. 0. _ ~ 0. 0* 0. 0#' 0. .2298-02 

•222E— 010* 0. 0# 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 

RESIDUAL MEAN IS .61450E-15 ' _ ... . ... - 

SIGMA SO OF RESIDUAL IS .79732E-04 


• 559E+09 FOR VARIABLE 


PARTIAL F 
PARTIAL F 
VARIABLES 
PARTIAL F 
PARTIAL F 


1 

IS 

IS 

1 

IS 

IS 


MAXIMUM F VALUE IS 
VARIABLES IN REGRESSION 1 
VALUE FOR VARIABLE 
VALUE FOR VARIABLE 
IN REGRESSION 1 
VALUE FOR VARIABLE 
VALUE FOP VARIABLE 
PERCENT VARIATION FX?LAINED IS 100. CO 
S TO • DEVIATION OF RESIDUALS IS .382E-06 
TOTAL F VALUE IS .59773E+1C 
NEW PARAMETER ESTIMATES AND STO. DEV. ARE 
,_.70oE+000. • 5CCE-010 . C. 

« SIC E-050 • .2C3E-C50. 0. 

RESIDUAL MEAN IS .fr9970E-15 
SIGMA SO OF RESIDUAL IS .58433E-11 


0 0 0 

•11QE+11 
•546E+Q9 
0 0 0 

• 118E+11 
.5A6E+09 


3 

0 


0. 

0. 


o.. 

0. 


0* 

0. 


•630E-04 


-4i 


2 
-si 

ss. 

fo trr 

- <0 tj “ 

c g 
> o 
£T pi 

— f 

■< 09 





RUN 1 EQUATION 1 NPTS 43 



THE 

V 

ALPHA 

0 

DEIS 



• 19 OC 6 +01 

• 3500 E + 02 

.5374E-01 

. 24202+00 

-.1732E+00 


.19506+01 

.35006+02 

.65376-01 

.2647E+00 

-• 1 73 2E + 00 


• 2900 E +01 

.35006+02 

.63506-01 

-.31736+00 

-.51756-01 

— 

.29506+01 
_ .A650E+01 

• 35006 + C2 
.35006+02 

.51706-01 

.35536-01 

-.35656+00 

•1953E+00 

— * 38 8 4E-01 
-.39846-01 


• 4 700 6 + 01 

. 3 5C0E + 02 

•4322E-01 

.19196+00 

-.39616-01 



•47506+01 

.35006+02 

.50476-01 

• 18 68 E + 00 

-.39176-01 



•4c Q0E+C1 

.350CE+C2 . 

. 5 722E-C1 

.17966+00 

39066-01 


.46506+01 

• 3 5 C0E + 02 * 

.63416-01 

♦17116+00 

-.3 8 f 46-01 

— 

.49OCF+01 
* 5 * 5C 6 +01 

.35006+02 

•350CE+02 

♦68976-01 

.64316-01 

.16066+00 
— • 1 0 2 1 E + 00 

-♦ 3 86 1 E-01 
-.36146-01 

' 

• 55 CG6 + 01 

. 3 5 00 E + 02 

.57606-01 

-.1290E+00 

-.35706-01 



• 55 5G 6 + 01 

• 3 5 00 E + 02 

.49816-01 

-.15616+00 

-.3 53 66-01 


•56006+01 

•35CCE+02 

.40666-01 

-.18336+00 

-.34686-01 


« 71 CO E +01 

• 35 CC6 + C2 

.35306-01 

.13096+00 

-.37826-01 



• 7 15 CE 01 

.35C0E+02 

.46076-01 

.12376+00 

33396-01 

U1 

.7200 6+01 

• 35 ( G6 + C2 

.44356-01 

•1153E+CO 

-.39066-01 

NJ 

.725uS+01 

• 35 CO E + 02 

• 4 81G 6 -0 1 

•1C56E+GO 

-.39506-01 


.73006+01 

• 3 5 COE + 02 

• 51 2 7E -01 

• 9 4 69 E-01 

-.3984E-01 


•735CE+C1 

.35006+02 

.53846-01 

•82566-01 

-.40066-01 

. 

.74005+01 

• 35 CO E + 02 

• 5577E-0 1 

.69296-01 

-.4039E-01 

_ 

•74505+01 

• 3500 £ + 02 

.57C26-01 

. .54866-01 

-• 402 86-01 



.75GC6+C1 

.3500E+C2 

.57576-01 

.39366-01 

— . 402 8 E-0 1 

i 

•75506+01 

• 35C0E + 02 

.57396-01 

.22896-01 

-.40396-01 



.76006+01 

.35C0E+02 

.56466-01 

.547^6-02 

-.40176-01 


•76505+11 

• 35 CO E + 02 

•5476E-01 

-.12756-01 

-.4CG6E-01 



• 7 700 E +C1 

• 3 5 00 E + 02 

.52306-01 

-.31646-01 

-.4028F-01 

- - 

• 77 5 0 6 + Q1 
. 78 00 E + 01 

.35006+02 

.35006+02 

.49066-01 
• 4 504 E-0 1 

-.51146-01 
— • 71 0 7E-01 

-.40066-01 

-.40176-01 


.78506+01 

.35006+02 

.40266-01 

-.91356-01 

-.40066-01 


.92O0E+01 

.35006+02 

•3555E-01 

* 8 79 3 E-C 1 

-.3434E-01 


•925CE+01 

.35006+02 

.38346-01 

• 774 2 E-0 1 

-. 3445E-01 


♦ 9300 E +01 

.35006+02 

.40526-01 

• 6 5 54 E-01 

— . 333 4 E-01 

— 

. 93 5 C £ +01 
.94006+01 

.35006+02 

.3500E+02 

•4208E-01 
. 4 298 E-0 1 

.52476-01 
_ .38336-01 

-.33006-01 
-* 32 e 3 E-01 


.94506+01 

•35C0E+02 

.43216-01 

•2328E-01 

-.33226-01 



.95C0E+01 

•35C06+C2 

.42726-01 

.7163E-C2 

-.32436-01 


.95506+01 

.35006+02 

.41516-01 

-.97916-02 

-.32326-01 


.9*006+01 

.35CCE+02 

.39566-01 

— • 2 7 39F -0 1 

-.32666-01 



.965CE+01 

•35C0E+02 

.36916-01 

-.45616-01 

-.32666-01 



.11006+02 

.35006+02 

.36406-01 

•556C6-C1 

— .3547E-01 

[ ■ 

. 11055+C2 

• 35 CC 6 + 02 

.37726-01 

•4323E-C1 

-.35136-01 

r 

• 11 1 C F +02 

.35006+02 

.36436-01 

.29846-01 

-.34' U “01 



APPENDIX 3 


This appendix contains the output generated by PROGRAM STEP for 
example 1 . 
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R T ’ - 0 o C* 

FCRT « •5E+01* 
XTRIKOP « 1, 

IPSKP - 10, 

L AT DP » 0, 

I AC E LOP - 0/ 
IFILOP * C, 

SEND _ _ 1 . 


TRIM VALUES 

ALPHT 8ETT AILT DELET OELKT 

• 2095CE+00 0# __0. -•83180E-01 0# 


T3 O 

0 2 

2 > 

33r“ 

O "0 

. c > 

> 6> 

m 

“3 _• 

< W 



SINPUT 


X PR E SOP 


1# Zt 3# 


•1374E+02, 

•996E+01# 


• 1 4 E + Ol* 


• 10 55 E + 04# 

• 1C272P+C1* 


•981E+01# 


• 23 57 E +04# 
•3051 E + 04* 


• 48 33 E *04* 


• 177E +03 # 
-.83186-01* 


•2095E+00* 


■ C. 0 p 


APPENDIX 2 


This appendix contains the NAMELIST/INPUT/ for example 1. 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


r~ 

! 

\r 




J 

r~ 


f" 



i 


a si 3 
Nani3d 

_*9efifi9/(G33SI)iV01d-k0UNVa s 

tillin' QKV*Q33SI-a33SI " " 

SV£9T*Q33SI*G33$I 

_ /I*ZZe/G33SI ViVQ _ 

(r)WOQKVd NQIlDMAi ~ X 


X 30 Vd ZZ*T* # IO •8Z/0T/28 8ZS+8**Ni3 X*,UG VZ/VZ. WOONV* NOI-lONGd 



SUBROUTINE PSET 74/74 OPT-X . ' FTN 4.8 + 528 82/10/28. 07.41.22 


X 


3 


xo 


c 

c 

c 


15 


4> 

o> 


20 


25 


30 


35 


C 

C 

C 


SUBROUTINE PSET 

COM NON /FLAGS/ X PSK P, NPTS , IOXH, JDIH, NNAX, I MIN, ICNT < 25 ) , IORO ( 25 ) 

♦ , IPPTS,XTRIH0P,ICAU 

COMMON /START/ X ( 25, 500 ) , XX SUN ( 25, 25 ) , XSUNI25 >, X3 AR ( 25 ) , S IGMA ( 25 ) 
- COMMON /ORDER/ IEO,N 

DIMENSION ITEMP(500),XTEM(5C0 j 
IRAN-0 

SET INDEX FOR SELECTED DATA PTS 
IPPTS-0 

DO 5 4-l,NPTS, I°$KP 

IPPTS-IPPTS+X 

4P«J 

IF ( IRAN .HQ. 1) JP-RANDOM< J)*NPTS 

I F ( J P .EO. C) JP*1 ' 

I T E M P { I PP T$ ) ■ JP 
5 CONTINUE 

DO 1C 1*1, N 

DO 20 Jol»XPP7S _ ~~ ” 

XTEMtJ )-X(I,ITEMP( J)) 

20 CONTINUE 

CO 30 J ■ 1 , IPPTS 

X(I, JJ-XTEMU) 

30 CONTINUE _ ^ _ 

10 CONTINUE 

DO AC K-1,625 
XXSUM(K)-0.0 
40 CONTINUE 

DO 50 Il-lt N 

DO 5 C I *1, N " ” 

DO 50 J-l, IPPTS 

XX5UM( I# 1 1 )*XXSUM( I# II ) + X( I# J )*X( II*«J ) 

50 CONTINUE ~ 

REDEFINE 4DIH FOR REDUCED R OF DATA PTS 



PAGE 


40 


JDIM-25+IPPTS 

RETURN 

END 




1 SUBROUTINE REDE F 

’ COMMON / FLAGS/ I PSKP* NPTS* I DI N# J D IH#NMAX> IMIN* ICNT < 25 > # IORD( 25 ) 

* * IPPT$*ITRI*OP#XCAlt 

' CONNON /ST APT/ X(25* 500)# XXSU* ( 25 > 25 > # X$U*< 25 > * XB AR ( 25 ) * $ IGrtA t 25 ) 

5 C OfiMON/ACP / A(Z5*525>>0(5C05#PPRT(25)*PRSNIN 

_ ' COMMON /ORDER/ I E Q * N 

c 

C INITIALIZE A,0,I0RD 

- c 

___ 10 NP1-N+1 . 

DC ICO 1*1* N 

DO ICO J«1*N . . _ 

ACI,J)*XXSUM(I,J) _ .... 

100 CONTINUE , .... . _ ... 

15 " DC no 1*1. N _ . 

K*0 ■ . _ 

DC 110 J *NP1* JDIH ... 

4 >__ ' K-K+l ‘ _ . 

Ln " " A(I*J)»X(I*K) _ 

20 _ 110 CONTINUE _ : 

DO 120 I *1* IPPTS _ 

0(1) -0.0 _ _ 

120 CONTINUE _ _ ; _ . ...• *. 

130 CONTINUE ... ... __ 

25 DO 140 I *1# N • 

IOR D ( I ) * I . .... . ^ 

™' 140 CONTINUE .. 

RETURN _ .... _ . ; 

END ; „ .. ....... .. .. 


. — J 
l 

- i 

* t 

. .. j 










DO 350 1 * 1 , N 

CO 350 J-NP1»JDIM 

*5 AU, J)-AP(I» J) _ 

... . . ..... 250 CONTINUE 

' C UPDATE 0 ' ‘ * 

C _ ^ ' ■ ■ - 

50 DO 400 I«l> IPPTS " ’ 

C(I)*0(I)+A(IDEF/N*I)*A(IDEF»N*I> 

_ 400 CONTINUE 

' 10 CONTINUE 

_ RETURN 

__55 END 



I ORIGINAL PAGE I? 
! OF POOR QUALITY 


SUBROUTINE UPDATE 


74/74 OPT" 1 


FTN 4*8>52S 


82/ 10/26 • 07 *41 *22 ^ PAGE 1 


1 


5 


10 


15 


fO" 


20 


23 


30 


35 


40 


SUBROUTINE UPDATE 

CCMMON/FL AGS / I P $KP,NPTS# IDXN, JDI M, NMAX* XMIN, ICNT « 25 ) , IORD ( 23) 
♦ , 1PPT$>ITRIN0P,ICAIUIACE10P,IFU0P 

CONMCN/AOP/ A< 25,525 >,Gt 500) > PPRTC25) , PRSNIN 
COMMON /ORDER/ IEO,N 
DIMENSION APC25,525) 

C 

C UPDATE A MATRIX <«APR) USED IN PRESS CALC 

C . 

NM1-N-1 $NP1*N+1 

CALL REDE F . . 

KK-0 _ 

DO 10 11*1, NM1 - .... 

IF(ICNT(II).E0*0> GO TO 10 _ 

KK* KK ♦ 1 • _ . . . 

IPEF-KK 

DO 210 Ll-1,NM1 . ... . 

IFCIGRDUU.EO.m K-IL 

210 CONTINUE 

CALI INTRCHG(IDEF,K) „ __ 

K-IDEF . . . ... 

DO 100 1*1, IDI M . _ _ 

DO ICO J-l, JOIN . . . . ... 

AP<I,J)*A(I,J) . ...... 

__ 100 CONTINUE 

DO 2 CO I*K # N .. _ . 

DO 200 J*K, JDIM . . . 

IF C I.EO.K) GO TO 225 . . _ 

APU9J)-A(IjJ)-(A(K,X)*ACK»jn/A(K*K) . 

CO TO 200 _ 

225 CONTINUE 

IF( J.EO.K) GO TO 250 _ 

AP( I# J )-A(I, J )/SORT< A(K,K> ) 

GO TO 2C0 .... . 

250 CONTINUE 

A P ( I , J ) "SORT ? A ( K,K ) ) _ 

200 CONTINUE 

DC 3C0 I-lfN . . 

DO 300 J ■ I, N ... 

AU» J)-AP(I#J) 

A<J,I)«A(I,J) 

300 CONTINUE 




- -I 


r 





SUBROUTINE PRESS 

COMMON /FLAGS/ I PSKP>NPT$> IDI K# JOIN* NMAX, IMIN, ICNT( 25 ) » IORDI25 ) 

* »IPPTS>ITRIMOP,ICALL»IACEL0P. tFILOP 

COMMCN/START/ Xl25>500)»XXSUM<25*25)>XSUK(25)>XBAR(25)jSIGMA(25) 
COMNON/AOP/ A{25>525)»0(500)»P?RT(2J)»PRSMIN 

COMMON /ORDER/ I EQ» N 

NM1-N-1 

IOEF-O 

IMIN-0 

PRSMIN-1.E06 
DO 10 I*1>NM1 

ICEF*IOEF»ICNT(I> . ... — 

10 CONTINUE _. .. . 

INITIALIZE A»a, IORO IF IOEF-O . 

IF ( IDEF .GT. 0) GO TO 130 . .. . • 

CALL REOEF ...... — 

130 CONTINUE . 

COMPUTE PRESS (-PRSMIN) FOR VARIABLES NOT TET IN MODEL 


00 200 I I-l » NM1 

IF ( ICNTIin.EQ.lJ GO TO 200 

00 210 LL*1»NM1 . . ... . . 

IFMOROCLD.EQ.II) K-LL 
210 CONTINUE 
PRS-C.O 

on 2E0 I-1»IP®TS 

PNt)M-A(K,K>*A<N»I*N)-A(K»N>*AtK#)UI> 

PONM-A(K,K)*(l.-Ot I))-( A(K,N»I)*A(K»N+l>r 

PRS-PRS+l PNUH*PNUM ) / ( PDNH*PDNH) 

250 CONTINUE 

VRI TE ( 6# 8C0 ) PR S# 1 1 

800 FORMAT (5X» C PRESS* *E12 .5> 5X» *FOR VAR IABLE*> 15) 
PPPT ( t I ) »PRS 

IF(PPSHIN .LT. PRS ) GO TO 200 
PPSMIN-PRS 
IM.IN-II 
200 CONTINUE 
PETUPN 
END 
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85/05/02# 13.03.28 


PAGt 


11 


227 

229 



116 

196 

203 

204 


88 


125 


148 




program stepspl 

73/74 

OP T * 1 


STATEMENT LABELS 

DEF LINE 

references 


7014 

402 


14 9 

146 


0 

410 

INACTIVE 

16 1 

160 


0 

.411 


153 

151 

152 

0 

418 

INACTIVE 

164 



7070 

419 


163 

159 

160 

0 

42 0 

INACTI Vt 

168 

167 


0 

421 

INACTIVE 

169 

168 


0 

426 


179 

17 7 

178 

7132 

427 


175 

172 


7135 

42 8 


176 

170 

171 

7157 

429 


180 

166 

167 

7041 

430 


155 

183 ... 


7163 

431 


184 

157 

15 8 

0 

439 


191 

190 


bln . 

443 


136 . 

134 


6762 

444 


139 

137 


7246 

448 


206 

202 


7250 

44 9 


20 7 

201 

205 

6714 

450 


12 1 

235 

239 

0 

451 


210 

209 


Q 

JibQ . 


219 

215 


0 

461 


218 

217 


7342 

470 


231 

220 


6352 ^ 

800 


58 

47 


0 

801 

INACTIVE 

50 

49 


6325 

802 


52 

49 


6333 803 


... 5 4 

Jt9 


6347 

80* 


56 

46 

51 

2 6*27 

805 


68 

57 


6424 806 

..... .... . 

_ 67 

JS0 - 

61 

0 

807 

INACTIVE 

60 

59 


6373 

800 


62 

59 


6410 

809 


.65 

59 


7732 

950 

FMT 

240 

139 


7741 

951 

FMT 

24 1 

155 


7746 

952 

FMT 

24 2 

16 2 


7755 

953 

FMT 

243 

169 


7762 

954 

FMT 

244 

19 5 


7770 

95 5 

FMT 

245 

197 


10013 

956 

FMT 

249 

213 

214 

7713 

95 7 

FMT 

200 

199 


7776 

956 

FMT 

Z*t 6 

212 



FTN 4.8052 
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PAGE 
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m 
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16 4 
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STATfc M E 
~10004 " 
10016 
7361 
0 

7361 

0 

10023 
7345 
7352 
~ 6 
7364 
6226 
Q — 

LOOPS 
‘ 6233 
6251 
6260 
6276 
6315 
6163 _ 
6450 
64 56 

64 67 

6500 
6505 
6515 

65 23 
6524 
6535 
6547 
6555 
6564 
6572 
6603 
6613 

—5625 

6635 

66 50 
66 64 
6671 
6703 


PROGRAM STtPSPL 

73/7^ 

□ P [ ■ l 



NT LABELS 

DEE LINE 

.^FE^ENCtS 


959 

eht 

247 

35 



99 6 

FMT 

250 

138 



.1000 


252 

40 


232 

15 OC 

INACTIVE 

25 3 




1999 


251 

135 


140 

2000 

INACTIVE 

135 

134 



2001 

FMT NO Kt FS 

25 4 




2111 


234 

135 



2112 


236 

233 



2500 

INACTIVE 

255 




2501 


257 

31 



2502 


27 

256 



2 503 

INACTIVE 

3 2 

31’ 




FTN 4.8+552 


35/05/02. 13.03.26 


PAGE 
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LABE L 

49 
30 

1000 

604 

606 

206 

51 

52 

53 

50 
50 
10 0 
100 
100 
101 
101 
201 
200 
202 
202 _ 

203 

204 

204 

205 
210 
301 


INDEX 

FROM- 

TO 

LENGTH 

PROPERTIES 



11 

28 

29 

58 

INSTACK 



I 

33 

34 

46 

INSTACK 



I 

35 

35 

4B 


EXT 

PEES 

L 

40 

25 2 

1 066B 


EXT 

REES 

I 

4 8 

5b 

3 4B 

OPT 



I 

58 

67 

4 3B 

OPT 



I 

74 

75 

3B 

INSTACK 



I 

76 

77 

2B 

INSTACK 



I 

78 

79 

3B 

INSTACK 



I 

30 

81 

36 

INSTACK 



I 

82 

85 

153 


NOT 

INNER 

IT 

84 

85 

26 

INSTaCK 



II 

Jb 

09 

2 38 


NUT 

INNER 

I 

87 

89 

20b 


NOT 

INNtR 

J 

6 8 

89 

3B . 

I NS T AC K.„ 



II 

90 

92 

1 46 


NOT 

INNER 

J 

91 

92 

36 

INSTACK 



I 

93 

96 

168 


NOT 

INNER 

J 

94 

95 

3b 

INSTACK 



I 

97 

100 

218 


NOT 

INNER 

J 

98 

100 

5B 

INSTACK 



T 

101 

102 

76 


EXT 

KEFS 

1 

103 

10 7 

238 


NOT 

INNER 


NUT INNER 


105 

108 

110 

114 


10 7 
109 
112 
115 


56 

36 

6B 

26 


1NSTACK 

INSTACK 


EXT PEES 


INSTACK 



PROGRAM STEPSPL 


FTN 4.8+552 


65/05/02 


73/74 CPT * i 


13.03.26 


LOOPS 

LAoEL 

index 

FROM- 

TU 

LENGTH 

PROPERTIES 







6774 

401 

II 

144 

150 

30B 


NOT 

INNER 






7010 

401 

JJ 

14 5 

150 

1 IB 

OPT 








7025 

411 

12 

151 

153 

14a 


NOT 

INNER 






7032 

411 

J 2 

15 Z 

153 

3u 

instack 








7052 

419 

11 

15 9 

163 

2 16 


EXT 

REFS 






7075 

429 

11 

166 

16 0 

6 5B 


EXT 

RtFS i 

NOT INNER 




7110 

426 

II 

170 

176 

326 


NOT 

INNER 






7126 

425 

JJ 

171 

176 

11b 

OPT 








7143 

426 

12 

177 

179 

14B 


NOT 

INNER 






7150 

426 

J 2 

175 

179 

3b 

INSTACK 








7174 

439 

II 

190 

191 

3B 

INSTACK 








7226 

449 

II 

201 

207 

2 5B 


EXT 

RtFS 






72 56 

451 

11 

209 

210 

46 

1 NS TACK 








7302 

460 

II 

215 

219 

2 26 


NOT 

INNER 






7313 

461 

JJ 

217 

218 

3b 

INSTACK 








COMMON BLOCKS 

LENGTH 

ML MB 

ERS - 

BIAS NAME ( LENGTH ) 









START 

37720 


0 

X 

(36000) 


36000 

XXSUM 

(1600) 

37600 

X SUM 

(40) 





3 76 40 

XBAR 

(40) 


37680 

SIGMA 

(40) 





AC DATA 

18 


0 

SAREA 

(1) 


1 

6 SPAN 

(1) 

2 

C BAR 

(1) 





3 

M 

U) 


4 

RHU 

(1) 

5 

G 

(1) 





6 

IX 

(l) 


7 

IY 

(1) 

8 

12 

(1) 





9 

1X2 

(1) 


1C 

OELET 

(1) 

11 

ALPHT 

(X) 





12 

BETT 

(1) 


13 

0 ELAT 

(1) 

14 

OELRT 

(1) 





15 

QT 

(1) 


16 

PT 

(1) 

17 

RT 

(1) 


FLAGS 

46 


0 

NPTS 

(1) 


1 

ICNT 

(40) 

41 

LATQP 

(1) 





42 

ITRIMOP 

m 


43 

ICALL 

(1) 

44 

1 ACE LOP (1) 





45 

IF1LUP 

(i ) 









ORDER 

2 


0 

IEQ 

(i) 


1 

N 

(1) 





„ ACCtL 

9000 


0 

AX 

(900) 


900 

AY 

(900) 

1800 

A2 

(900) 





27 00 

PDOT 

(900) 


3600 

QDOT 

(900) 

4500 

ROOT 

(900) 





54 00 

VEL 

(900 ) 


6300 

P 

(900 

7200 

Q 

(900) 





01 00 

R 

(900 ) 









KNOT” 

23 


0 

XNOT 

(23) 








EQU I V 

CLASSES 

LENGTH 

ME MB 

LRS - 

BIAS NAM 

t ( L t Nb Trt ) 









A 

6241 


0 

RR 

( 62t 1 ) 








STATISTICS 













FROGK AM LENGTH 

4 

54 57 B 

19247 









BUFFER LENGTH 



5b 70 B 

3000 









CM 

LA Bt IF 0 COMMON LENGTH 

13 

33 31 B 

46509 
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3DM 1N0D 009 

0C9'Z0S < (T)V1VT 19M Si) )(I TOOfc 
TO 08 '9 666 ( ( T) 303) 31 
(hdn'wmp )vivo) (Davjs zus 

OOT-maT* X S * ( X > d 1 1 T*NT1 XOHH E- ( X )dOH * 0Z-H9NT T»OI 

TONIiNOD TOS 
ZOS'6666 ( ( T ) 303 ) 31 

aaH'(HDN'I*I'(I)SlINni)'(HDN‘T«I'(I TSVWVNJ'HDN'OI ( T)OV 3d IOS 

T-Sidr*SidN ir * 03 • { Sidr-sidN) ) ji 

z r*z»sidr 
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96 Ui 09 (T *i 9* 11 VDI ) dl 
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(006) d'(0C6)0'(006)d' (0 06 IDA ‘IT 
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iHd IV '131 30 '7 XI '7I'AJ 'Xl'9'OHX'W'D'fl'S /V IV OD V/ NUWWUD 
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(006)3130'(006)X13(1'(006) VllO N0ISN3WT0 
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(Z)XI NUTSN9WT0 
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(X 'A'i 'Si)l 3SV1V0 DNIinoXflnS 
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r* 

SE 

OE 

sz 

oz 

ST 

01 

c 
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3 N 1 1 1 10 d q OS 



SUBROUTINE DATASET 7J/7<t OPT-1 


FTN 4.B+552 o5/05/0.2» 13»03#28 


45 


50 


53 ' 


60 


65 


-~4 

oo 


70 


75 


80 


601 CONTINUt 

CALL GETRANt 1R,IN>2»RN»Y1» Y2> 

1N-1N+1 

T ( I ) ■ 0ATA<1)-TS 
VL L ( I ) * DA TA (2 ) 
atTA(I)- DATAti8)-BETT 
ALPHU t* 0 AT A( 19) 

P( I ) * DATA t5 J 
Q(l)-DATAt6)+.003*RN 
Rill* DATA (7 ) 

AX t I )* DATA 11 0) 

AY ( 1 )■ DAT At 11 ) 

AZ(i)» DAT A( 12 ) 

DC L A ( I ) * D AT At 13 ) -DEL A T 
DTLRtl)- DATAt 15)-DELRT 
DELL (1 >*DATA (14J-DELET 

PDOT ( I ) ■ DA TA 1 16) SQDOTt l> *DATA< 20) 1R0 OT 1 1 ) * D AT A ( 1 7 ) 


15 CONTINUE 

IMIACELOP.EO.O) GO TO 46 

CALL SECDERt 3# 3>T> P#POO» PT EM * NPT S» PO » P3> .5>. 5» WRK1*WRK2> WRK3» 
1WPK4 jWRK5» WRK6 ) . „„„ 

CALL SECDER(3»3»T»Q#QDD»QTEH»NPTS*P0>P3#.5*.5»WRK1*WRK2»WRK3» 
1WKK4»WRK5»WRK6) 

CALL $ EC UE Rl 3» 3t T > R» ROD* RT EM »NP TS» P0> P 3, , 5, . 5> WKK1> WRK2» WR K3 , 
1WRK4,WRK5> WRK6 ) 

DO 45 I* 1> NP TS 

PDOTU ) • OE RSPtTUlxTiPxNPTS* POD,PTEH ) 

QDQTtl )*DERSPlTt I)»T*Q*NPTS#QDO>QTEM) 
RDOT(I)-DLRSPtT(I),T,R,NPTS,RDD,RTEH) 

45 CONTINUE 

46 CONTINUE 

IF(LATOP.NE.l) GO TU 803 
DO 800 I *1 *N PT S 
Xll,I)*BtTAtI> 

X(2,I)-P (1 )* B/ (2.*VEL( I) ) 

Xt3/1)*R(1 )* 6/ 12. * VE L ( I)) 

X(4,I)«DLLAtI) 

X(5> l)-DELR(I) 

X ( 6> I ) “ALPHt 1) *X 1 1 # I ) 

Xt 7, I) «ALPHt I) *Xt 2»I ) 

Xl8,I)»Xt3>l)*ALPHtl) 

Xt 9> 1) *UEL At I) *ALPHt I) 
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( I'2)X»U 'T?)X 1 1 6) 10NX *3 9* (I JHdIV ) dl 
It UONX-II )H«nv*(I ‘02 >X l|6>10NX*39‘ IIIMdlV) dl 
« I‘2JX«II ‘6T)X <<mONX*39* mHdIVUT 
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1 1 ‘7 1 X» I I ‘/IX I (?) 10NX *39* (I IHdIV) dl 
l7T10NX-{T)Hd1VV(r r 9TX II7)I0NX*39* IT IHdIV HI 
U‘2)X-(I*S)X (IT) 10NX *39* (I )Hd3V) dl 
I 1 1 10 NX-1 I)Hd3V»( I‘6)X (IT) 10NX ’3 9* (I IHdIV) dl 

•o-li'Tinx s oh 
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SU BP DiJT I N E 


133 


135 


1*0 


1*5 


150 


1 55 


160 


165 


DATASET 73/7* 0PT*1 


FIN *.8+552 85/05/02. 13.03.28 PAGE * 


IE ( A L P H ( 1) .GE. XNOT (10) ) X ( 22 » I ) « AL PH ( I )-XNOT (10) 

IF (ALP Hi I) .GE.XNOT U01 ) X < 23 , I > * X ( 2, 1 ) 

IF (ALPHt I ) .GE. XNOT (11) ) X{ 2* , I ) * AL PH ( I )- XNOT ( 1 1 ) 

IF ( ALPH< 1 ) .GE. XNOT (11) ) X(25,I)-X(2, I) 

IF I ALPHt 1) .Gt. XNOT (12) > Xt 26, 1 ) » ALPH < I l-XNQT ( 12 > 

IF ( ALP HI I).GL.XN0T(12t) X( 27, 1 1 *X< 2, I) 

IF (AtPH( I) .GE. XNUT(13) ) X ( 28 , I ) » AL PH 1 1 )-XNOT (1 3) 

IF I ALP H< I) .GE.XNOT (13) ) X( 29, I ) *X(2, I) 

IF (ALP H( I) .GE. XNOT (1*) ) X( 30, I ) * AL PH ( I I-XNOT (1*) 

IF ( ALPH( I) .GE. XNOT(IA) ) X ( 31 » I ) *X ( 2, 1 1 

IF ( ALPH( 1) .GE.XNOT (15) ) X( 32, 1 ) * AL PH ( 1 )- XNOT < 1 5 > 

IF ( ALP H( I) .GE.XNOT (15)) X( 33, I > * X{ 2, I ) 

IF (ALPHt I) .GE. XNOTtib) ) X( 3* > I ) • AL PH ( l )-XNOT (16) 

IF < ALPHt I) .GE. XNOT (16) ) X( 35 > 1 ) -X( 2, 1 ) 

IF (ALPHt I) .GE. XNOT (17) ) X( 36, I ) * AL PH (I >- XNOT ( 1 7 > 

IF (ALPHt I) .GE. XNOT (17) ) X( 37 , I ) »X l 2, 1 ) 

IF (ALPHt 1) .GE.XNOT (7)) X( 38, I ) «X( 3, I ) 

801 IF (ALPHt I) .GE.XN0TU1J ) X( 39, I ) «X( 3, I ) 

80* CONTINUE 

PRINT 96*, 10,IE0,NPTS 
IF (1CALL.GT.1) GO TO 999 
IF (LATOP .EQ.O) GO TO 50 
PPINT 968 

96 8 FOKMAT (IX, 7X,* TIME ♦, 11 X, *)/♦, 12 X, *8ET A* , 1 IX , * P* , 1 *X , * R*, 1 IX , 
1*L>ELA*,11X,+ 0ELR*) 

96* FORMAT (1H1,3X, *RUN*» 15, 5 X, +EQUATIQN *, 12 , 5X, *NPTS* , I 5, / / / ) 

PRINT 96 2, (Ttl ),VEL(I),BETA( I),P(I),R(I),DELA( I),DELR(I),I*1,NPTS) 
962 FORMAT (7 (2X, El 2.* ) ) 

CALL INFOPLT (1,NPTS, T, 1, DE LA , 1, 0 . , 0. ,0., 0. ,0., 1, 1HT, 

1*, *HOE LA, 2 2, 7. ,5., .75, .75) 

CALL INFOPLT (1,NPTS> T, 1, DE LR , 1 , 0 . , 0 . ,0., C.,0., 1, 1HT, 

1*, *HL)E LR,2 2, 7. ,5. , .75, .75) 

GO TO 999 
50 CONTINUE 
PRINT 998 

998 FORMAT (lX,7X,*TIML*»llX»*tf*»llX,*ALPHA*» llX,*0*,l*X,*0tLt*) 

PRINT 99 7, (T(i),VEL(I)»ALPri( I),Q(I),OL lE ( 1 ), l*i,NPT$) 

997 FORMAT (5 (2X,E12.*> ) 

CALL INFOP LT (i,NPTS, T, 1, Ot Lfc, 1, 0 ., 0. ,0 . , 0. ,0 ., 1, 1HT, 

1*, *HUF LE,2 2, 7. ,5. ,.75, .75) 

GO TO 999 
) 99 8 PRINT 9999 
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FUNCTION uERSP (XX#X>Y,N>P#H) 

FUNC T I UN DERSP 

UCR$P IS USED TO UiiTAiS THE FIRST DERIVATIVE OF SPLINE 
CURVE E ITTED DATA 

USAGE - 

X * DERSP (XX# X*Y#N*P#H) 

NOTE - IF XX LESS THAN X i 1 1 THEN DERSP * UXtl)/DY 

IF XX GREATER THAN X(N) THEN DERSP - DX(N)/OY 

WHE RE — 

XX f NO fc PEN Dt NT VARIABLE FOR WHICH INTERPOL ATEO ! 
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N-DIMENSI ONEu VECTOR OF 1ND EPENOt NT POINTS 
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$U 8RCIU T 1 Nfc SECDERC Ll,L2#X, Y, P,H,N, P0,P3, XKI, XK2, A, B,C,D, GaMMA, 

1 BETA) 

SUBROUTINE SECDER 

SPCOtR IS USED WITH 1 UNCTION SPLINE TO PkEPQkM A SPLINE 
INTERPOLATION. IT IS USED TO GENERATE P A NO H. 

USAGE - 

CALL SECDER(L1, L2, X, Y, P, H,N,P0,P3, XKI, XK2, A, 8, C,D, GAMMA, BETA) 
WHERE - 

Li,L2 DETERMINE THE END CONDITIONS AT XU) AND X(N) TO BE 
USED. (SEE BELOW! 

X N-D 1MEN SHINED VECTOR OF INDEPENT POINTS 

Y N-D IME NSIONED VECTOR OF DEPENDtNT POINTS 

P N-D 1MEN SI ON ED VECTOR TO BE RETURNED 

H l N-l) -DIMENSION ED VECTOR TO BE Kfc TURNED 

N NUMBER OF DATA POINTS 

StCUND DERIVATIVES ARt GIVEN AT THE END POINTS 
XK1 NOT USED 

XK2 NOT USED 

IF LI *1 THEN 

PO SECOND DERIVATIVE AT X(1),Y(1) 

IF L2-1 THEN 

P3 SECOND DERIVATIVE AT X(N),Y(N) 

FIRST DERIVATIVES ARE GIVEN AT THE END POINT 
XKi NOT USED 

XK2 NOT USED 

IF LI *2 THEN 

PO FIRST DERIVATIVE AT XU), YU) 

IF L 2* 2 THEN 

P 3 FIRST DERIVATIVE AT X(N),Y(N) 

NO INFORMATION ABOUT THE CURVE IS KNOWN 

PO NOT USED 

P 3 NOT USED 

IF 11*3 THEN 

XKI P 1 * ( 3 , 0 ) * XKl*P*»(3,i), XKi GREATER THAN 0 

IF 12*3 THEN 

XK2 P»M3,N| * XK2* P • M 3, N-l ) , XK2 GREATER THAN 0 

A N -U IME N Si ON ED WORK VECTOR 

■i N-D IME N Si ON ED WORK VECTOR 

C N-D IME N SI ON ED WORK VtCTUR 
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D N-D 1ME N SI ON ED WORK VECTOR 

BETA N-D IM t N Si ON fcD WORK VECTOR 

GAMMA N-D IMENSlONfD WORK VICTOR 

SUBROUT XNES CALL - 

none 

DIME NS ION X( N> f Y ( N )f A ( H) t 3 (N ) t C ( N) *D <N GAMM A( N ) p BET A( N ) #H (N )> P ( N ) 
K-N-l 

DO 1 J *1 # K 

1 HU)«X(J+1 )-X( J) 

DO 2 J*2>K 

A( J ) - 

B ( J ) -2.* (H(J ) + H ( J- 1 ) )/H( J) 

C( J) * 1. . 

2 D( J ) = 6./H(J)*((YU*l)“Yl J) )/HU)-(Y(J)-Y(J-i)l7HlJ-U.) 

IF (L1.EQ.2 ) GO TO 20 

IF (LI. F Q .3 ) GO TO 10 
8(11-1. 

C C 1 1 *0 • 

D(1)-P0 
GO TO 30 

10 B ( 1 ) *1 • 

C( 1) — XK1 
D ( 1 ) «D • 

GO TO 30 

2 0 BU1-HUJ/3. . . 

C( l)«H(l)/6. 

D ( 1 ) •( Y l 2 ) -Y (1 ) ) / H (1 1 - PO 
30 IP (L2. EQ.2 ) GO TO 21 
IF (L2.EQ.3 ) GOTO 11 
Ai N ) -0 • 

B ( N ) *1 • 

D ( N ) -P 3 
60 TO 40 

11 MN)->XK2 
‘ 8(N)«1 . 

D ( N ) -0 * 

GO TO 40 _ . 

2l“ A(N j«H(K j/6. 
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•o 

TO-3 8Z.P2 * 

00+39992* 

20+3C0SE* 

20+302 TT * - 

•o 

T0-35*88* 

00 + 16552 * 

Z0+3C0SE* 

20 + 35XXX * 

•o 

TO-3 0iO8* 

00 + 3T6 5? * 

Z0+300SE * 

20 + 30T TT * 

•o 

T 0-3 T6 88 * 

00+18292* 

20+30098* 

20+350 TT * 

•o 

TO-3 CfcZ.2 * 

00 + 3 PT 92 * 

20+30098* 

20 + 300 TT * 

•o 

TO-36588 * 

00 + 3 Z6 62 * 

20+30098* 

20+3 96 OT * 

•o 

XO-XSe<>E • 

OO+3T062 * 

20+30098* 

20+3060T * 

*0 

To-ieeoe* 

00 +3 69 62 • 

20+30098* 

20 + 3590X * 

*0 

TO-3 6968 # 

00+3 Lb 62 * 

20+30098* 

20 + 3 Ofl OT * 

•o 

TO-359T** 

00+36262 * 

20+30098* 

20 + 35Z. OT * 

•o 

TO-3 68 08 * 

00 +3 IT 6Z * 

20+30058* 

20+30Z0I* 

•o 

TO-38258* 

00+38682* 

20+30098* 

20+3 590T* 

•o 

TO-3 2898 * 

00+35Z82 * 

20+30098* 

20 + 309 OT * 

•o 

T 0 -3 8288 * 

00 +3 £9 £2 * 

20+30058* 

20 + 359 OT * 

•o 

t o-Sfwtr* 

00+30* e~z * 

20+30098* 

20 + 305 OT*' 

•c 

TO-3 9Z62 * 

00+382 €2 ‘ 

20+30098* 

20+356 OT * 

•o 

io-3oeoe* 

00 + 3 80 82 * 

20+30098* 

2O + 30Vgn 

•o 

T0-30628 * 

00+38622 * 

20+30098* 

20 + 3 58 OT * 

•o 

TO-3<^9T * 

00 +36Z 22 * 

20+30098* 

20 + 308 OT * 

•o 

T 0-3 906T * 

00+3 99 22 * 

20 + 3 CO 98 * 

20 + 3 52 OT * 

•o 

T0-3T69T • 

00 + 3*9 22 * 

20+3C098* 

20+302 OT * 

•o 

T0-3886T • 

00+3 2* 22 * 

20+30098 * 

20 + 3CTOT * 

•o 

T0-3995T * 

00+328 22 * 

20+10098 * 

20 + 3 OT OT * 

•o 

2O-3Z90Z * 

00 +3 22 22 * 

20+^0098 * 

20 + 1 50 OT * “ 

•o 

20-3 9Z. T? * 

00 + 3*1 22 * 

20+31*098 * 

20+3C0CT • 

•o 

20-3Z598 • 

00+3 90 22 * 

20+30098* 

10+30566* 

•o 

20-? °062 *- 

00 +3 TO 22 • 

20+30098* 

TO +3 00 66* 

•o 

80-3686 1 

00 + 3 96 T 2 • 

20+ 3 CO 98 * 

10+30596* 

•o 

20-15650 •- 

00+3 86 T2 • 

20+300 98 * 

10 + 30096* 

*0 

■ "20-TP9t>6 *- 

■“ 00+3 26 12 * 

20+ j 0098 * 

TO + 305 Z 6 * " 

•o 

TO-3 6^ 2T *- 

0O4326T7 • 

20+30098* 

T0+300Z6 * 

•o 

T0-12E8T *- 

00+3*612* 

20 + 30098 * 

10 + 3^996* 


_ MAXIMUM F VALUE IS .3A9E+03 FUR VARIABLE 3 

VARIABLES IN REGRESSION OOIOOUOOO 

VARIABLES IN REGRESSION OOOoOOOOO 

PER(LENT_VAR1ATI0N EXPLAINtD IS D9.A9 

STD. DEVIATION OF RESIDUALS IS ,70bE-0l 

TOTAL F VALUE IS .3AR01E+03 

NEW PARAMETER ESTIMATES AND STD, DEV, ARE 

0. 0. -.231E+010. 0. 0. 0. 

0. 0* 0* 0* 0, 0, 0, 

0* 0. 0. 0, 0. 0, Q, 

0. 0. #124 E +000 • 0. 0. 0* 

0 • 0* _ . . 0 • 0 • 0 , . .0 • 0 • 

o* 6* o • o • o • o* o* 

RESIDUAL MEAN IS -• 4631 3E-1 5 

SIGMA. .SQ OF RESIDUALS IS .11819E+Q1 


MAXIMUM F VALUE IS •691E+03 FOR VAR I A8 Lt 1 

VARIABLES IN REGRESSION 101000300 

VARIABLES IN REGRESSION _00Q000000 

""PARTIAL F VALUE FOR VAR I AbL t 1 IS •6B8E+03 

PARTIAL F VALUE FOR VARIABLE 3 IS .109E+04 

_VARlABL|S_iN REGRESSION 101000000 

VARIABLES IN REGRESSION OOOOOOOOO 

PARTIAL F VALUE FOR VARIABLE 1 IS . 698E+03 

PARTIAL F VALU E FO R VARIABLE 3 IS . 109E+04 

PERCENT VARIATION EXPLAINED IS 89,65 
t- STD, DEVIATION OF RESIDUALS IS .358E-Q1 
S TOTAL E val ue IS ,102:5E*09 

NEW PARAMETER ESTIMATES AND STD. DEV. ARE 

-.lOSE-FOlO. -.209E + 010. 0. 0. 0. 

0. O'* o, 0, 0, 0, 0. 

0, 0* 0* 0* 0 • 0* 0* 

^WE^piQ*.' , 62 ^ E-01 0, 0. 0. 0. 

o, o, 6, o* o • o»- o • 

0, 0, 0, 0, 0. 0, 0. 

_ RESIDUAL MEAN IS ,797351-13 
'SIGMA"5g OF RESIDUALS IS .30185E+00 


MAXIMUM F VALUE IS .575E+G2 FOR VARIABLE 2 
VARIABLES IN REGRESSION 11 10CO 0 0 0 

VARIABLES IN REGRESSION 0 C C 0 0 ) ) C U 


ooo oooooo 


OODOOOoOOOOOOOO 
O 0 0 0 0 0 


0, 0, 0, 0, 0, 0, 0, 0, 

0, 0, 0, 

0, 0, Ot 0, 0. 0, 0, - • 723E-0 

0, 0, 0. 0, 0, 0* 

0, 

0, 0, 0, 0, 0# 


oooooooooooo 

0 0 0 


oooooooooooo 

ooo 


0 , 0 . 0 . 0 , 0 , 0 , 

c. 

0* 0, 0, 0, 0, • 246E+0i 

0 , 0 . 0 . 0 , 0 , 0 , 

0 . 

0, 0, 0, 0, 0, 


o 0 0 0 0 0 0 0 0 0 0 0 0 0 0 

o g 6 o o o 


ORIGINAL PAGE IS 
OF POOR QUALITY 



PARTIAL 

F 

VAL UE 

FOR 

VAPIABLE 

i 

IS 

. A96E+03 



PARTIAL 

F 

VALUE 

F UR 

VAR 1ABLE 


IS 

. &73E+02 



PARTIAL 

F 

VALUE 

FOR 

VARIABLE 

3 

IS 

.12<t_E+04 



VARIABLES 

IN REGRESSION 1 

1 

1 

0 0 0 0 

6 

0 

VARIA8L tS 

IN R E GR F S S I ON 0 


0 

0 0 0 J 

0 

0 

PARTIAL 

F 

VALUE 

FOR 

VARIABLE 

I 

IS 

• <t96 E +03 



PARTIAL 

F 

VALUE 

FUR 

VARIABLE 


IS 

• ^7 3fc +u2 



PARTIAL 

F 

VALUE 

FGR 

variable 

3 

IS 

. 12 4 1 +0*» 




PERCENT VARIATION EXPLAINED IS 91.68 _ 

STD. DEVIATION OF RESIDUALS IS .321E-01 

TOTAL F VALUE IS .B6327F+03 

NEW PARAMETER ESTIMATES AND STD. DEV. ARE. _ 

-•908E+00— .683 c +oi— .236E+010. 0. 0. 0. 

0. 0. 0. 0. 0. 0. 0. 

0. _ 0. 0. 0.. Q.. . _ 0, 0». 

.357E-01 .715E+00 ,56<*F-010. 0. 0. 0. 

0. . 0. 0. 0. Q. 0. 0. 

0. 0. 0. C. 0. 0. 0. 

RESIDUAL MEAN IS .68295E-13 

SIGMA SO Of RFSIDUALS IS .2A27U+09 - . 


MAXIMUM ( Value IS . 103L+03 FOk VARIABLE 2a 
VARIABLES IN REGRESSION 111000000 
VARIABLES IN_REGR ES S I ON 0 0 0 1 0..._Q ... 0 .... 0. 0 

PARTIAL F VALUE FOR VARIABLE 1 IS .362E+03 

PARTIAL F VALUE FOR VARIABLE 2 IS . 569E+02 

PARTlAL_t -VALUE JOR ..VARIABLE.. 3 IS -.119E+09 - 

'PARTIAL F VALUE FOR VARIABLE 20 IS .103E+03 

^ VARIABLES IN REGRESSION 11 1 0 0 0 J 0 0 

o VARIABLES IN REGRESSION 0 0 0 1 0 fl . Q_ 0 0 

00 PARTIAL F VALUE FOP VARIABLE 1 IS .362E+03 

PARTIAL F VALUE FGR VARIABLE 2 IS .56AE+02 

PARTIAL F VALUE FOR VARIABLE 3 IS .U9t+09 

PARTIAL F VALUE FOR VARIABLE 28 IS .103E+03 

PERCENT VARIATION EXPLAINED IS 9A.22 
STD. DEVIATION 0F_ RESIDUALS IS .268E-01 . . 

TOTAL' F VALUE IS ^95 390E+D3 

NEW PARAMETER ESTIMATES AND STD. DEV. ARE 

-.730E+00-.572E*01-.211E+010. 0. 0, 0. 

0 , 0. 0. 0. C. 0. 0. 

0. 0. -. 1 AOE+OIO. 0. 0. 0. 

0. 

0. 

0. 


.29BE- 01 .597E+00 .A71E-010. 0. 0. 

0 . 0 . 0 . 0 . 0 . 0 . 
0. 0. .105 E + OCO. 0. 0. 

IKLTiHIAI Ml- /.Al fS .*><,1791-13 


o o o o o o 


ooooouuoooooooo 

0 C O 0 0 0 


0. 

• 209 E + GG 

0. 


00000000000000 
0 0 0 0 0 


oooooooooooooo 

o a . , o _ .... o 


0. 0. 0* 0* 0. 0* o. 

0. 0, 0. 

o. o.. 0. o«_. .. a. o • o. 

0 , 0 . 0 * 0 . 0 . 0 . 

0. 0. 

o. c. 0 • G • 0. 0. 


0# 0. 0§ 0 • 0* 0* 0* 
0.0. 

()• 0. 0. 0« 0* 0« #l73E + 00 

c. 0. 0. 0. 0 . 0 . 0 . 0 . 

0. 0. u. 

0 . 0 . 0 . 0 . 0 • 0 « 0 . 



SIGMA SO OF PfSIDUALS IS .IbBSdc+OD 

maximum ^ value is .ssal+oe for variable is 
variables IN REGRESSION lllOOJOOO 

VARIABLES IN REGRESSION OOOIOOOOO 


PARTIAL 

F 

VALUE 

FOP 

variable 

1 

IS 

. 2156*03 

PARTIAL 

F 

VALUE 

FOR 

VARIABLE 

2 

1$ 

. 119E *03 

PARTIAL 

F 

VALUE 

FOR 

VARIABLE 

3 

IS 

. 1036*0* 

PARTIAL 

F 

VALUE 

FUR 

VARIABLE 

15 

IS 

. 5526*02 

PARTIAL 

F 

VALUE 

FOR 

VARIABLE 

2 8 

IS 

. 179E *03 

VARIABLES 

IN REGRESSION 1 

1 

1 

0 0 l 

VARIABLES 

IN REGRESSION 0 

0 

0 

l 0 < 

PARTIAL 

F 

VALUE 

FOR 

VARIABLE 

l 

IS 

4 21 5 6 *03 

furniL 

r 

VALUE 

FOR 

VARIABLE 

2 

IS 

•1196*03 

PARTIAL 

F 

VALUE 

FOR 

VARIABLE 

3 

IS 

.1036*0* 

PARTIAL 

F 

VALUE 

FOR 

VARIABLE 

15 

IS 

•552E+02 

PARTIAL 

F 

VALUE 

FOR 

VAR IABL E 

28 

IS 

• 1.796 *03 


PERCENT VARIATION EXPLAINED IS 95.33 

STD. DEVIATION OF RESIOUALS IS .2A2E-01 

TOTAL F VALUE IS .950B9E+03 

NEW PARAMETER ESTIMATES AND STD. DEV. ARE 

-.583L+00-. 128E+02-.193E+010. 0. 0. 0. 

0 . 0. 0. 0. 0. 0. 0. 

0. 0. “.186E+010. 0. 0. 0. 

•269E-6I .538E+00 .A25E-010. 0. 0. 0. 

0 » 0 . 0 . 0 . 0 . 0 . 0 . 

0 . 0 . . 9496 - 010 . 0 . 0 . 0 . 

RESIDUAL MEAN IS .43231E-13 
5 SIGMA SO OF RESIDUALS IS .1363OE+0O 



MAXIMUM F VALUE IS 

. 369E+02 FOR VARIABLE 

31 


VARIABLES 

IN REGRESSION 1 

1 

1 

0 0 0 

0 

0 

0 

VARIABLES 

IN REGRFSSION 0 

0 

0 

10 0 

1 

0 

0 

PARTIAL F 

VALUE FOP VARlAbLE 

1 

IS 

. 262E+03 




tartial f 

value for variable 

2 

IS 

. 132L+U3 




PARTIAL F 

VALUE FOR VARIABLE 

3 

IS 

• 109E *0* 




PARTIAL F 

VALUE FOR VARIABLE 

15 

IS 

• 71 7E *02 




PARTIAL F 

V AL UP F Ok VARIABLE 

2 8 

IS 

. 11 66 *03 




PARTIAL F 

VALUE FOR VARIahI E 

8 1 

IS 

. 36 7 E *02 




VARIABLES 

IN REGRESSION 1 

x 

X 

0 0 0 

0 

0 

0 

VARIABLES 

IN REGRESSION 0 

0 

0 

10 0 

1 

0 

0 

PARTIAL F 

VALUE FOP VARIABLE 

1 

IS 

. 262 1 *03 




PARTIAL F 

V AL U - FOR VARIABLl 

2 

I s 

• 13 2 £ *u 3 




P ART 1 al F 

VALUE 1-Ht. variable 

S 

IS 

. 109 F *0* 





o o o 



o. 0. 0. 0. 0. 0. 0. .107E+C 

0 * 0 • 0 • 

0 . 0 . 0 . 0 . 0 * 0 * 0 . 41386*0 

04 0 . 0 . 0 . 0 . 0 . . 6506*0 

0. 04 

0« 0. 0. 0. 0. 0* 


00000100 0 000 0 0 0 
0 0 0 C 0 0 


0 0 0 0 0 0 0 


0 C 0 

0 C 0 


0 0 10 0 
0 0 0 


ORIGINAL PAGE IS 
OF POOR QUALITY 



•o 

00 43T85 * *0 

00 + 30TT * *0 

20430ET* *0 


•o *0 *0 

•o *0 *0 

•o *0 *0 

•o *0 *0 


•0 *0 *0 

•0 *0 *0 

•0 *0 *0 

•o • 0 *0 

•o *0 *0 

•o *0 *0 


0 0 0 0 0 0 
OOOOTOOOOTOOOOO 


0 0 0 0 0 o 

OOOOTO OOOTOOOOO 


*0 *0 *0 *0 ’0 *0 *0 

•o *0 *0 

00 + 1 509 • *0 *0 *0 *0 *0 *0 *0 

00+3ZVT* *0 *0 *0 *0 *0 *0 *0 

•0 *0 *0 

2049 * TT « *0 *0 *0 • 0 *0 *0 *0 


# 01 04 3601 * *0 *010-3**** *<J *o 

•o *0T0-12T** *0 *0 *0. • o 

•0 *0 *0 * 0 TO-3 09E ' 0043Id** TO-3T*2* 

• 02 Q+ 320T "" '*0 *000+3 126*- r 0 *0 

•0 *000436X5*- *0 *0 *0 * 0 

•jO *0*0 *0T0+38T2*-2C+32*T *-00+ 352**- 

3d V *A 30 * (US QNV SiivwlIsT miWVdVd «3N 

E0+3VUS8* $1 3m v A 3 3Y1QX 

. . .. . T0-3 9T?* si - 

0E*96 $1 aiNIV1dX3 NQUVldYA lN33»3d 




204U5* * 


SI 

TE 

3 *19 V I dV A dUJ 303VA 

3 ivnavd 



2049SE2* 


SI 

92 

3 18V I dV A dO 3 3I33VA 

3 3V t x dV d 



2043502 * 


ST 

02 

3iflviavA aod 3nivX 

3 IVUdVd 



2049686 * 


SI 

5 I 

3 18V i dv a yu 3 amvA 




*043021* 


SI 

E 

3 38V I dVA dO 3 3fnvA 

3 IVUdVd 



e043£9T* 


SI 

2 

nsvIavAaud 3nTvA 

dlvilavd - 



2043**9 * 


SI 

T 

339VIHVA dO 3 3fnVA 

3 IVIldVd 

0 

T 

0 0 T 


0 

0 

0 NOl$S3d93d NT 

$319VTdVA 

0 

0 

0 0 0 


T 

T 

T NO I SS3 d9 3d NT 

S319VIHVA 



204 32. 5* * 


SI 

TE 

319VIHVA MLH 3f11VA 

3 3VT IdV d 



HO* 35 £2 * 


SI 

02 

3iaviavA aud smvA 

J 



2041502 * 


SI 

02 

3 39V I dV A d03 303VA 

3 IVIidVd 



H0+36E6* 


SI 

5 T 

3 is? t av a aod amvA 

3 IVTldVd 



V0+30HT ♦ 


SI 

f 

3 39V I dV A dO 3 303VA 

3 IVIldVd 



E0+3E9T * 


SI 

7 

3 39V I dV A d 0 3 SHIVA 

3 IVIidVd 



HO* 3* V9 * 


SI 

T 

3 39 V I dV A dU 3 3D 3V A 

3 iVUdVd 

0 

T 

0 0 T 


0 

0 

0 NO I S $ 3 d9 3d N I 

S3 38 V IdV A 

0 

0 

0 0 0 


T 

T 

T NO I SS3d93d NI 

S3 38 VI dV A 

02 

319VIHVA HOd ? 0* 35 02 * 

SI 3D3VA 3 WflWTXVW 





oo+Booiii* si sivnaisia do 

OS VW9I£ 







£T-19T0<?V SI NV3H 

3VndlS3d 

OT043ETT 

• • 

0 



*010-362.8* *0 

*0 

0 


*0 

0 



*0 *0 *0 

*0 

0 


•o 

0 



• 0 TO- 3 56£ * 0043TC5* 

T0-3052* 

0X043**6 

• _ • 

0 



•0T049E5T*- *0 

*0 

0 


*0 

0 



*0 *6 *0 

*0 

0 


*0 

0 



•0TO432I2*-2C43S2T*- 

0043X09*- 


3dV *A3(J *OIS ON V SUVWI1S3 d313WVdVd KIN 
F 0+ 3 T T 026 * S I If! IV A 3 fVlOX 
TO-3522* SI STVOOISiM 30 NOTi VTA 30 *{US 
2-6*5 6 SI U3N I V3 d X -i NOIiVldVA iNjlOdld 
204 3Z9E * SI Te 31PVI4VA dOJ 3niVA‘ 3 IVIldVd 

€ 04 3 8 T I * S I 02 33HV)dVA dO 3 303VA 3 3VIIdVd 

204 3Z.T/. * ST qj 338VIHVA dU J JOTVA 3 OVTiMVd 



RESIDUAL Mr AN IS , 330U2;.-13 

SIGMA SO HE ^FSI DUALS IS .10PCtlE+OO 



_ j 

MAXIMUM F ' 

VALUE IS 

t 15 3E +0 2 

! HUP VAKIABLc 

21 

VARIABLES 

IN REGBFSSION 1 

1 

1 

0 0 0 

0 

0 

VARIABLES 

IN P- 

GRES SI ON 0 

0 

0 

10 0 

1 

0 

PARTIAL 

F 

VALUE 

FOR 

VARIABLE 

1 

IS 

• 701 E +0£ 



PARTIAL 

F 

VALUE 

FOR 

Variable 

2 

IS 

. 175 E +03 



PARTIAL 

P 

VALUE 

FOR 

variable 

3 

IS 

.1Q2E+04 



PARTIAL 

P 

VALUE, 

FOR 

variable 

15 

IS 

. 355E+02 



PARTIAL 

F 

VALUE 

f OR 

VAR I ABLE 

20 

IS 

. I68E ♦02 



PARTIAL 

F 

VALUE 

f OR 

VARIABLE 

21 

IS 

•152E+02 



PARTIAL 

F 

VALUE 

FOR 

VARIABLE 

28 

IS 

. 336E+0* 



PARTIAL 

P 

VALUE 

FOR 

Variable 

31 

IS 

.557E+02 



VARIABLES 

IN REGRESSION 1 

1 

1 

0 0 0 

0 

0 

VARIABLES 

IN REGRESSION 0 

0 

0 

1 0 .0 . 

. 1 

0 

PARTIAL 

F 

VAL UF 

FOR 

variable 

1 

IS 

•701E+02 



PARTIAL 

F 

VALUE 

FOR 

VARIABLE 

2 

IS 

. 175E+03 



PARTIAL 

F 

VALUE 

FOR 

VARIABLE 

3 

IS 

. 102E+04 



PARTIAL 

P 

VALUE 

FOR 

VARIABLE 

15 

IS 

• 35 5 £ +02 



PARTIAL 

P 

VALUE 

FOR 

VARIABLE 

20 

IS 

. 168E+02 



PARTIAL 

F 

VALUE 

FOR 

variable 

21 

IS 

•152E+02 



PARTIAL 

F" 

VALUE 

FOR 

variable 

2 8 

IS 

. 336E+Q2 



PARTIAL 

F 

VALUE 

FOR 

VARIABLE 

31 

IS 

• 557E+02 




PERCENT VARIATION EXPLAINED IS 06. 62 
STD. DEVIATION OF RESIDUALS IS .2I0E-O1 
TOTAL F VALUE IS .79856E+03 
_NfcW P AS AM cj LB .ESTIMATES AND STD. DEV. ARE 


-.A30E+00-.IA3E+D2-.200E+01D. G. 0. 0. 

i- 0. 0. 0. 0. -.A6IE+0O .60DE+010. 

£ Q. Oa_. -.110E+01D. 0. . -• U.OE +0 20 • 

.233E-0I .A67F+00 .369E-010. 0. 0. 0. 

0. 0. 0. 0. .'.OOE-Ol , 6A 2£ tOOO . 

0. 0. . 3 2DE-01 0 • 0. . lODE *010 • 


RESIDUAL MEAN IS .33A70E-13 

SIGMA $y OF RESIDUALS IS .10I3bE+0C 


■ • - ■ ■ • — 

M A K 1 HUN 

F 

VALUE IS 

. 10 OE +02 

FOR VARIABLE 

29 


VARIABLES 

IN PEGPFSSrON 1 

1 

1 

0 C 0 

0 

0 

0 

VARIABLES 

IN REGRESSION 0 

J 

0 

1 1 0 

1 

0 

0 

PARTIAL F 

VALUE 

F OR 

VAR IABLE 

1 

IS 

• 709E+02 




PARTIAL F 

VALUE 

FOR 

VAR I ABLE 

2 

IS 

• 182E+03 




PARTIAL F 

V ALIN: 

FOR 

VAR I ABE E 

3 

IS 

. 10 IE +04 




PARTlAt f 

value 

^ OP 

VAP IAULE 

15 

is 

. 374E+02 





o o o 


0 C 0 0010 0001 1000 

U 0 0 0 0 0 


0 C 0 0 01000011000 

oooooo 


0. 0. 0. 0. 0. 0* .947E+0 

0 . 0 . 

o* 0, 0 ... 0* 0. 0. •llOE+O 

04 0« 0* 0. 0. Ot 0* • 564 E ♦ 0 

0 * 0 * 0 * 

0 . 0 . 0 * 0 « 0 . 0 . 0 . 


OOOOOiOOOO 110 0 0 

0 0 0 0 0 0 


ORIGINAL PAGE IS 
OF POOR QUALITY 



PARTIAL 

f 

VALUE 

ERR 

VARIABLE 

2 u 

IS 

• 207 1 +02 

PARTIAL 

F 

VALUE 

FUR 

VARIABLE 

2 x 

IS 

. 19GE+u2 

PARTIAL 

F 

VALUE 

FOR 

VAR I ABLE 

28 

IS 

. 250E+02 

PARTIAL 

F 

VALUE 

F OR 

VARIABLE 

29 

IS 

. 100L+02 

PARTIAL 

F 

VALUE 

FOR 

VARIABLE 

31 

IS 

. 771E+01 

VARIABLES 

IN REGRESSION 1 

1 

1 

0 0 0 

VARIABLES 

IN REGRESSION 0 

0 

0 

1 1 0 

PARTIAL 

F 

VALUE 

FOR 

VARIABLE 

1 

IS 

, 709E+02 

PARTIAL 

F 

VALUE 

FOR 

VAR I ABLE 

2 

IS 

. 132E+03 

PARTIAL 

F 

VALUE 

F OR 

VARIABLE 

3 

IS 

• 101E+0A 

PARTIAL 

F 

VALUE 

FOR 

VAK IABLE 

1 5 

IS 

. 37AE+02 

PARTIAL 

F 

VALUE 

FOR 

VARIABLE 

2 0 

IS 

. 207E+02 

PARTIAL 

F 

VALUE 

FOR 

VAR I ABLE 

2 1 

IS 

.190E+Q2 

PARTIAL 

F 

VALUE 

FOR 

VAR I ABLE 

28 

IS 

. 250E+02 

PARTIAL 

F 

VALUE 

FOR 

VARIABLE 

29 

IS 

• 100E+02 

TaETIAL t 

'value' 

FOR 

VARIABLE 

31 

IS 

• 77 IE ♦01 


PERCENT VARIATION EXPLAINED IS 96.67 

STD. DEVIATION OF RESIDUALS IS .206E-01 

TOTAL F VALUE IS .73873E+D3 
NEW PARAMETER ESTIMATES AND STD. Dfc V. ARE 
-.425E+00-.143E+02-.216E+010. 0. 0. 0. 

0. 0. 0. 0. -.506E+00 .6646+010. 

0. 0. -.96IE+00-.694E+010. -.5976+010. 

72296-01 .458E+0(T.362E-010. C. 0. 0. 

0. 0. 0, 0. .3926-01 .630E+000. 

0, 0. .804 £-01 .904E+000. .103E+010. 

RESIDUAL MEAN IS" .32614E-13 
SIGMA SO OF RESIDUALS IS .971416-01 


- MAXIMUM 

F .] 

VALUE IS 

,610t +01 FOR VARIABLE 

13 


~ VARIABLES 

IN REGRESSION 1 

1 

i 

0 0 u 

0 

0 

0 

VARIABLES 

in Regression o 

0 

0 

110 

1 

0 

0 

PARTIAL F 

VALUE 

FOR 

VARIABLE 

1 

IS 

• 53<»E+02 




PARTIAL F 

VALUE 

FOR 

VARIABLE 

2 

IS 

• 16 9 1 +0 3 




PARTIAL F 

VALUE 

FOR 

VAR IABLE 

3 

IS 

• 103E+0A 




PARTIAL f 

VALUE 

FOR 

VAR IABLE 

13 

IS 

.607E+01 




paRtTSIT f 

VALUE' 

FOR 

VAR IABLE 

15 

IS 

. 697E+Q1 




PARTIAL F 

VALUE 

FOR 

VARIABLE 

20 

IS 

. 250E+02 




PARTIAL F 

VALUE 

FOR 

VARIABLE 

21 

IS 

• 195E +02 




' PARTIAL' f 

V AL De 

FOR 

VARIABLE 

28 

IS 

. 29^E +02 




PARTIAL F 

value 

FOR 

VARIABLE 

29 

IS 

. 10 3 E +02 




PARTIAL F 

VALUE 

FOR 

variable 

31 

IS 

•785E+01 




variables 

In Regression i 

1 

1 . 

0 0 J 

0 

0 

0 

VARIABLES 

IN REGRFSSION 0 

0 

0 

i 1 * 

1 

0 

0 

PARTIAL F 

VALUE 

FOR 

VARIABLE 

1 

IS 

. 534E+02 




PARTIAL F 

VALUE 

FftP 

VAR IABLE 

? 

IS 

. 169 F +03 





o o o o o o 


COOGO IOOQQIIOOQ 
0 0 0 0 0 0 


0. 0. 0. 0. 0. .952E+01 

0. 

0. 0. 0. 0. 0. •109E+0Q 

o. 0* 0. 0. 0. • 553E + 00 

0 • 

0 . 0 . 0 . 0 • 0 . 


00310100001 1000 

0 0 0 0 0 0 


000101000011000 

oooooo 



PARTIAL 

F 

valuf 

F n* 

VARIABLE 

3 

IS 

.103E+0<7 

PARTIAL 

F 

VALUF 

FOR 

VARIABLE 

13 

IS 

. 60 7 E ♦Ci 

PARTIAL 

F 

VALUE 

F OR 

VARIABLE 

15 

IS 

• 697 E >01 

partial 

F 

VALUE 

F OP 

VARIABLE 

20 

IS 

.250E+02 

partial 

F 

VALUE 

FOR 

VARIABLE 

21 

IS 

• 195E *02 

PARTIAL 

F 

VALUE 

FOR 

VARIABLE 

26 

IS 

.2** E*02 

PARTIAL 

F 

value 

FOR 

VARIABLE 

29 

IS 

• 10 3 E +02 

PARTIAL 

F 

VALUE 

FOR 

VARIABLE 

31 

IS 

.785E+01 

PERCENT 

VARIATION EXPLAINED 

IS < 

?6. 76 



STO. DEVIATION Of RESIDUALS IS .20AE-01 
TOTAL F VALUE IS .68020E+03 

NEH_ PARANETtR ESTIMATES AND STO, DEV, ARE... ... ... 
-.384E«-dO-.158E«-D2-.2I6E+Olb. 0. 0. 0. 

0. 0. 0. 0, - ,5 61E + 00 , 665E+010. 

0* 0. -,939E+00-.698E'*'0i0., _ . 5966 ±010. 

• 227E-01 • R5 3E +00 .358E-010. 0. 0. 0. 

0, 0, 0, 0, _ .3B8E-01 .623E+000. 

6V" 0. .7956-03 .894E+000. .102E*010. 

RESIDUAL MEAN IS .30R3U-I3 

SIGMA SO OF RESIDUALS IS .9A620E-O1 


MAXIMUM F VALUE IS .396E+01 FOR VARIABLE 30 


<jO-- 


o o o 



o. C. 0. 0. .530E+C10. , 57AE*C 

0 * 0 * 

o. c. 0. _ 0. 0« 0. *1QIE*C 

o. 0. 0, 0. 0* • 5256+000* .5*7E+C 

0 . 0 . 0 . 

0. 0. 0. 0* 0. 0. o. 
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APPENDIX 7 


This appendix presents three optional forms for the model structure 
determination basis m SUBROUTINE DATASET. One longitudinal option and two 
lateral options are given. 

The longitudinal option given here is an example of the second-order spline 
which gives a smoother representation than the first or zeroth-order spline. 

Is usedTor tSTr P I* f °f * second -° rder «Pline in C When this basis 

°m ^nation, it provides a second-order splice for * 
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LONGITUDINAL EQUATIONS: SECOND OPTION 


an example for C z 


where 


°z ■ C z <0) 6 - , - S - 0 + c z U> <lc/2V + 

e q 


17 

C (a) = c (0) + C + T D (a - a ) + 
2 2 z a i - 1 V i 


C z «*) 
q 


+ C z <* + C 

qa 


z 2 

q® 


17 

+ I D (a - a.) 2 
i = 1 qa 2 . 1 + 


C (a) 
z 6e 



D. (a - a) 0 
6e i 1 + 


C (a) 6 
z 6e e 
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LONGITUDINAL EQUATIONS: SECOND OPTION 


ORIGINAL PAGE IS 
OF POOR QUALITY 


_bbj906_ W.NPTS 

x(i#ii-u>H(i) 

)*C/(2*VEL(I ) ) *0(1.) .. 
>13,1 f«0ELE ( I ) 


~xU»i )»ALPH(I ) —XKNOT (1) 


00 <507 III**. 3? 

907 X ( I II * 1 1 * 0 « 

IF(ALPH(I).Gt.XKHOT(l)) 

X ( 5. I ) *X< 2. I >* ALPHl I ) 

~IF ( ALPHt I ) . GE. XKNOT ( 2) ) X ( 6, I) • A l PH t I ) - XKNOT t 2 ) 

X ( 7, I) • X ( 2, I ) * A L P H (I) ** 2 — . :■ ■ - Y - ■ . . 

IFIALPHI I ). GE.XKNOT (31 ) __x_t 0 » I 1 *A LPH 1 1 )-XKSuT ( 3 

IFTaLPHI I l.GE.XKNOTU) ) 

IF ( ALPHt 1 ) < GE « XKHOT ( A ) ) 

IF (ALPH1I ).GE.XKN0T(3) ) 

•" ~_IFtALPH(n.6E.XKN0Tt5>> 


X [ u 9 \ (ML mi } — A a » U » A .» ( 

X ( 9» I ) * l ALPHl Il-XKNCT (I) 1 ♦ * 2 * X 1 2 » I ) 

X( 10» I > * ALPH ( I 1-XKNOT ( A) . ... 

Xtllr I ! "(ALPHt I ) - X K N 0 T (3! ) * * 2 ♦ X ( 2 » I ) 
IF (ALPHt I) .GE.XKNOT t 5> I X ( 12» I ) ■ ALPH l I J -X* HOT t 3 J 

IF I ALPHl 1 ) .GE. XKNOT ( 5J ) * | ^ ^ : i * K KO f T 7 7 ' ’ ’ ‘ ? ‘ ‘ ? * 

X(17»I)"lALPH(I>-XKHCTC6n*«?*x(2,n 
X ( 18» I). "ALPHt I)-X«.S0T(6 I 
X( 19, I ) • ( AtPHt I )-_>*NCJ t?n_**2*_x 12.11 
X ( 2 GViV- alp H t'n-X M* DT ( 9 1 

Xt 21* I ) ■ < ALPHt I)-X<SOTl r J 1*«?*x ( 2»I > 
X( 22/ I) • ALPHt I ) -XKNOT (101 
x'(23* I >~* (ALPHl I)-*XN0T(9 1 )*»2*X (2. 1 1 


I F ( A L P H ( I ) . G E • XKN 0 T ( 7 ) ) 
IFlALPHtl 1.GE.XKNQT (6) I 
IFIALPHI I) .GE.XKNOT (81 ) 

I F ( ALPH ( I) . GE . XKHOT m L 
'YfTaLPH(I j.GE. XKHOT (91) 
IFlALPHin.GE.XKNOTlSn . 
“IF (ALPH ( I) .GE.XKNOT 110)) 


i I :S«:;sSoI i »i i 

i f utpirriiTcfi XKHOT no) » x 1 2 5> i i * ( al ph j j irMNOi.t 10 
IF (ALPH ( I ). GE. XKNOT 1X2 1 )._* *26*1) "ALPH ( I ) -XKNOT ( 1 2 


TTfilPHf T \ GF. XKHOT tl2 } J XlcO# 1 I ■ A L r n i l i — a a n u i i 

- ? irKlf tl'li XKNOT 1 11 > r~X(27»X)"(ALP H 1 1 )-XKNOT ( U ) » ** 2*X t 2, I ) 

IfStfSrt ‘HU: ^ 

IF ( ALPH / 1 ) 1 ® *2*X 1 2 > 1 1 

x(32.n-ALPH(n-xKNom?» 

X(33*n*(ALPH(I1-XKN0T(lA)1**7*X(2.1) 

Xt3A»lL*ALPHII)-XKN0Tll6) 

Xl33»n*(ALPM(I)-X<NOT(15))**2*X(2.I) 

X(36,n*ALPH(I)-XKN0T(17) 

X 1 37 * I ) • t ALPHl I l-XKNOT 1 17 n ** 2 *x < 2 » 1 1 . 
X< 38 »I)-X( 3 >I) 

X t 39 » 1 1 "XI 3> 1 1 


906 


IF (ALPHl 1 ) . GE.XKNOT t 15) ) 
IFlALPHtn.GE.XKNOTmn 
IF ( ALPHl I ). GE.XKNOT (161 » 
|F(ALPH(I).GE.XKNOTtlJ) ) 
IF (ALPHl I ) .GE.XKNOT <171 ) 
IF jALPHlI 1. GE.XKNOT tl7) ) 
*1 Ft ALPHt I). GE.XKNOT (7) ) : 
IFIALPHI II. GE. XKNOT ( 131 ) 
CONTINUE 


/*> n ▼ r\ a a a 
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Next, we give 
simple version 
in B should be 


an example of a first option for the lateral equations. This 
incorporating zeroth-order splines in a and first-order splines 
used for the first approximations to the lateral coefficients* 
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LATERAL EQUATIONS: FIRST OPTION 


An example for Cn 


C n = C n (a> 6) p = r = 0 + C n ( a) P b / 2V + C n <a) rb/2V 

6=5=0 p r 

a r 


where 


+ C (a) 6 + C (a) 6 


n 6a 


a n 


6r 


C (a, 3) 
n 


y 

= c n + C 0 + l B . ( 3 - 0 ) 

0 1 i = 7 1 i + 


+Ca+ l C (a -a ) 0 

n a i - 1 "0 1 + 


C (a) = C + l C (a - a.) ( 


n n . “ , n 

p p i = i p. 


C (a) = C + l C (a - a ) ( 
n ^ n , n i 

r r i = 1 r - 


C n < a) ' C „ + £ V «■ - V 

6a 5a i = 1 6a . 


c (a) = C + l c (a - a )° 

n 6r n 6r i = 1 n 6r 1 + 

i 


and 


<6 - 6 t ) + 


+ ( 


0 for 0 < 6. 


3 - B i for B > 6. 

3 + 0 for B _< B 
i i 
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OF POOR QUALITY 


3flNI J.N03 L 08 
(I)Km><l'»CIX _ 

?6)iON»X*(t )fi3««ii 7 9C)X <t6)i.6NXX-*dl*(l)?138Hl 
<6)10NXX><im3t-U*«C).X ((6)iONXX*39’( 1)9139)31 

t9>10NXX+CnT13«-(I‘A£>X ltmONXX-'3V(im3liHl 

(mOWM-CimjS-tl'ZEJx < (B)10NXX*i9‘(n»13M)3r 

U )iONXX+( I)ti38-n <9t >X ( U JiUNXX-*3 1* C 1 )f 138 » 3X 

(i>10NXX-( 1) V139»< I*9£)X t U)10NXX‘39* (X 1*138)31 

U'SlX-CI'iElX <<9>10NXX*33M IlHdl*) 31 
Cl** IX* t I**£)X C C9 110NXX* 33* < IlHdT*) 3l 
U«£)X«CI*£E>X U9)10NXX , 39*t IJMdlT ) at 1 

' { PZIX-C I'ZE) X C (9110NXX* 39* C IlHdl* >3l 
C I*I)X«CI«X£>X « ( 9 | ION XX *33* (I )Hd1*)3l 
C I*S)X-CI*0£1X CCSllONXX'jgMllHdmdl 
<I**)X»II*6Z)X CCilXOHXX'ag'mHiTVHI 

ci ‘EJx-d'ezrt (CijiONxx'jg'mHdiyjdx 

I I*Z)X»U*1Z)X I ( 5 ) 10NXX* 39* (I)Hd1*>3l 
1 1 *t ) X •( I*9Z) X < Ci)10NjlX*33*CI>Hd1*Hl 
< I •i )X*CI*SZ)X CC*)10NXX*33*CllHdT*)3I 
C I**)X«CI*VZ)X CC<>)iDNXX*39*CI)Hd'l*)3l 
(|*£)X»CI*EZ)X CC*110NXX*33*mHd1*)3l 

CT*Z)X«C1*ZZ)X <C*U0NXX*39*CI)Hd1*)3l 
CI*I)X-CIMZ)X ( CyllONXX'dS'CI JHd1*)3l 
CI*<)X-CI*OZ)X C C£ )iUMXX*39‘( l)Hd1*)3l 
(I *8 ) X* ( I *6 1 ) X ((£)lONXX»39MIIHd1»> Jl 

<i*e>x»u*«iix ((euoMxx^g^ciiHdiyiiX 
( l*2)X»U*it)X ( U )10NMX*33*(1 IMdlf )3I 
ci*iix-ci*9tix c t euoNxx* j-j*mHdi«)di 

CI*S»X-4I*ST»X (CZ)iONXX*39*<l)Hdl*)3i 

ci**ixici**t>x cczjxoNdx'dgMnHdivcdi 

C J *E IX • C I *ET 1 X C C 2 ) ION XX *39*(l)Hd1*)3I 
C 1*ZIX-(I«ZT)X C tZ)10NXX*39*tI)HdlX)9l 

ci'nx-ci'uix ccziioNxx*39*mHdiy»3i 

U*t IX-CI'OTIX CCTU0NXX*39MIIHd1»)3l 
CI«V1X-CI*6)X C C t)10NXX*39*C l)Hd 1* 1 31 
1 1 *£) X ■< I *8) X < miONXX*39MI)Hd1V)3l 
U'Z)X»CI*Z)X (< t) ION XX * 39 MIlHd IX) 31 
C1*T)X»CI*9)X ( ClliONXX* 39* (I IHdl* )dl 

♦O-U'IIDX 06tf 
6E*9-III 668 DU 
C 1 ) Slid* ( I *4 ) X 
C I)*13U-C 1'*>A 
t< I) 13A**Z J/8«(Ild»C I*E)X 
( ( I) 13A«*< 1 / X* C 1 1 d* C I # Z » X 
CH»i3«-CX*T»X 
SldN*I*I 108 DO 

NOIldC IS8I3 iSH0Il»n03 1»*31*1 0 
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If the first option given for the lateral equations indicates a need for a 
two-degree spline in (a, 3), the following lateral option can be used. 
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LATERAL EQUATIONS: SECOND OPTION 


and example for Cn 

4 

C n (Q ’ 8) " c o + C 6 + l (A + A. 6) (a - a.) 0 

i = 1 i i 1 + 

7 4 7 

+ l 3 (S - 3 ) + l l D (0 - 3.) (a - a )° 

j = 6 °j J l i - 1 j = 6 ^ J I 1 + 

Note: for the analysis it was assumed that A q. =0, i= 1, 2, 3, 4. 

This assumption was confirmed by the later analysis using partitioned data. 


C (a) - c + £ C (a - a.*) 


*n0 


"p "P i = 1 V 


C n (o) = c „ + l C (a - a )° 


n r °r i = 1 V 


C (a) = C + l c (a - a )° 
n 6a n 6a i = 1 n 6 a 1 + 


C (a) = c + l . c (a - a )° 
"6r n 6 r i = i " 5ri i + 
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Ui -d 

Sa 3 

a. O’ 
.j o£ 

< o 

59 
o 51 


/ 


( 9) 10NXXM 1) VI 38«( I '>£) X . * 

UZI10N)(X‘« , mWH*0N»'(9liDN»I-'3r(I )1l3Q)il 
(9)10N>lx-U )*13U*CI **£> X * 

( ( Z )1UNV* *39“ { ])Hd1V QNV* (0 ) 1UNX* • 30 ° ( I ) *138 > 31 
U) iO*>U* < I ) Y138-C !.'££> » , * 

H!liU'l9> , 59 , ll)l.nrQN? , tmiNM-‘5T(llVlatiHi 

( l ) iomxx - ( i ) 7i nt • 1 1 t ) * * 

( ‘v) 1 1 1 v i <b-( i * zu x * 

( (l ) 1 ' • N V »*.'•) * l I ) V rll V* U r<v * { 6 ) 1 [> * » - * i l’t I) f i ’ d > -1 » 
(9)1JN>H-1 I ) V X jM « t i‘tl )A ♦ 

{ < U iUN'U • 3i> * M bid IV ONV I 9 )10N*1 • j-j" II>»13 « ) j 1 
tillJNm*U)UJ#-(I'UIX C (iUl<mx-* 3V (I >*138)Jl 
miUfldX-<I>*l?8M I Mt>X C lL )iUM*X*3 J* 1 i)*138> d I 
IVIlUNxxMI >* 138 -CI‘ 0 £>X t(<mUNx*-*3 l*(l)»l3«»dl 
C V ) 10NXX-1 1)*138*(I'0£)X (<9>1UMXX* J9MUV13UI3I 

{ -1*51 X«( I ‘6Z>X ( ( 5 )10NXX ' 39 * ( l)Hd1*):»I 

< I'*>X«C I««Z>X C (4 uONXX‘39*(l>H<i1*)iI 

<i'£)x-(i‘zz>x t csuosxx'aoMDHdiyjii 

( I 'Z> X •( I '9Z )x (<;>10HXX*39‘(l>Hd1*>3I 

C I '5 )X •( I ‘SZ >X ( t 6 UONXX* 39* 1 UHdTX) dl 
U'V)X«(I < VZ)X ( (9)10NXX*39* ( IlHdl*) 31 

( I '£ > X “( I # £Z > X ( (6)10NXX*39*(l>Hd1f)dI 

(I'Z JX«(I'ZZ)X 1 ( V UONXX* 33* ( 1 )rtd IT >31 
<I‘T)X«(IMZ>X l(V>10NXX*39*t I)Hd1*)dl 

(I‘t)X»CI'0Z>X C C£ > ION XX *39* t l>H«n*> 31 

n<V>X*II'6t>x‘ (~(£>10NXX*39 , Cl>Hd1*>3l 
<I‘E)X-CI'«l>X ( (£)10MXX*39*(I)Mdl*>il 

<I‘Z>X-tI‘it>X (U>A0rt»X*39*( I>Hd/)*>3l 

<I«X >X»U'9T>X ( C E HOHXX * 39 * ( 1 >Hd1*> i I 

(I'S>X»U'SX>X <(Z>10NXX*39MI>Hd1»>iI 

C I'6>X-t l«9I>X ( ( Z ) 1DNXXJ 39*11 >Hd1* > 31 

(I'e>X-(I'Ei>X C <Z)iONXX*39 , CI>Hd'»*>3i 

<I'Z>X-CI'Z1>X (CZI10NXX*39*(I)Hd'l*>3l 

(■IMIX*U‘UIX t <Z)10NXX*39* C I)Hd1*)3l 

( I ‘4 )X-( I ‘OX >X C miUNXX'd-J* C ISHd'lYMI 

(I'*F>X-CI'b)X ( ( T > 10HXX* 39*C I)Hd'l*)3I 

<I'£)X»tI'a>X (miUNXX*3V*C I ) Hd 1* > 3 1 

(I<Z)X*( l*X.)X { (1 1 10NXX *39*CX>Hd"l*)dI 
(I'l)X*CI«9)X ( miONXX*39*tI)Hd1V)31 

•Q-CI'IIDX 506 

' - ' • 6£ *VI1I 506 do 

11) 8130* (I'S >X 
(I>Y13Q*(I < Y>X 

• * T n > 1 3 A* ‘ z > / a « (I > a* n r t j x " 

CCI>13A* , <>/8*CI)d"CX , Z>X 

(1 ) V130-C 1 il>x 

Sid N '1*1 »0o 00 

NOIldO 0NO33S »SN011*no3 1*«31Y1 J 
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OF POOR QUALITY 







3DMI1N05 >06 

(~i ) lONXX ♦(I)V138»(I < 6E))( * 

t mi.ONXX’aO' t DHdU'aNV* U)iGNXX-*3V (UV139M1 
uuoNxx-mvi38-u'6E>x . 

UV>10NXX , 39MI)Hd'>**0Ny*(iJ10N)IX , 39*m*139>3l . .. 

(9)10NXXMI)*13B»U'8E>X * - 

( (>)1UN»X*39M nHd1V*0N¥* miONXX-»31MI >Ti38Hl 

t9)lONXX-(I)V138»CI'0£)X * ...» 

( ( VJlONXX'aO* ( I >Hdl**UN*‘ ( 9 > 10NXX • 39 * < I )*139>31 

(i)10HMX + mtl38»( l‘i£)X . , 

Ue>10NXX*39*U )H«nt*ON**(i>iONXX-*3VCl)*l39131 
tillONXX- (I)T138*I J *A E ) X * 

( (£)10NMX*39* ( 1 >H<n»*0MX* U UONXX* 39* « 1>*139 > 3.1 .. 

C9>16nXX + (I>»139-U‘^£)X m * 

t (£>iOHXX*39* ( 1 JHdlVQN** ( 9U0NXX-* 3 V U >*139>3l 


|9UONXX-(I>Vi39-«I'9£>X ....... . . , * , ... * 

( (£>iOHXX*3'j* ( I JHdlV'ONXMoiiONXx'ayM I)>13fl> 31 
U)10NXX + m*139»tI'«E>X * 

( (Z )XOHXX *39 ' ( I ) Hd 1 V *0N* * UU0NXX-*3 1* U 1*139)31. 
1l)iOHXX-£H*13«-ll*Sl)X . „ , . * 

t (d)10N>U*39* ( I IHdlV’ONV* ( iUUNXX ‘39* t 1 )*i.3*> dl 
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TABLE 


a 


<N 
Q Cl 


CO 

ca 


oa 


m 

ca 


CM 

a 

CM 

ca 


CM 

a 

ca 


a 

cx 


0 

u 


a 

ca 


a 

cx 


a 

M 


CM 

ca 


CM 

ca 

CM 

C 5 


cr 

a 


CO 

a 
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TABLE 2 


TRUE VALUE ESTIMATED ESTIMATED 

VALUE STANDARD ERROR 


c 

m 

o 

.105 

.101 


C m 

-.400 

-.384 

(. 023 ) 

a 




C m 

- 15.0 

- 15.8 

( 0 . 4 ) 

q 




Se 

- 2.00 

- 2.16 

(. 04 ) 

m ( A “)+, 13 0 

-.600 

-.561 

(. 039 ) 

m(Aa) + ^ I?,, 

0.00 

-.939 

(. 080 ) 

^(AcO^o 

- 1.00 

0.00 


m q(Aa)° +>9 o 

0.00 

5.30 

(. 52 ) 

Cm q(Aa)° +> i 0O 

+ 10.0 

5.74 

(. 55 ) 

c 

m / A \0 

q ( Aa) +}13<> 

+ 10.0 

+ 6.65 

(. 62 ) 

Cm q(Aa)° +170 

0.00 

- 6.98 

(. 89 ) 

Cm /A \0 

q ( Aa) +>180 

- 10.0 

- 5.96 

( 1 . 0 ) 
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a, deg 
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Figure 4.- Calcomp plotter output for STEPSPL in example 2. 
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TIME 

(b) Predicted and estimated pitching moment coefficient 
residual sequence, and autocorrelation function 
for one variable model. 


Figure 4.- Continued. 
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licted and estimated pitching moment coefficients, 
ssidual sequence, and autocorrelation function 
for three variable model. 
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Figure 4.- Continued 
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(g) Predicted and estimated pitching moment coefficients, 
residual sequence, and autocorrelation function 
for six variable model. 


Figure 4 .- Continued. 
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4 .- Continued. 
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(1) predicted and estimated pitching moment coefficients 
residual sequence, and autocorrelation function 
for eight variable model. 


Figure 4.- Continued. 
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(k) Predicted and estimated pitching moment coefficients 
residual sequence, and autocorrelation function 
for ten variable model. 


Figure 4 .- Concluded. 
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